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ABSTRACT 


This thesis discusses the implementation of a discrete 
control systems analysis option, ORACLS, for the CONTROLS 
analysis package. The option encompasses a collection of 
FORTRAN programs that facilitate the design and analysis of 
linear, multiple input/multiple output digital feedback 
control systems. All programs are interactive and operate 
in the VM/CMS environment of the IBM 3033. 


The main focus of this thesis was the development of the 
FORTRAN programs ORACLSX and TRANFUNC. They serve as the 
foundation of the ORACLS option, complimenting as well as 
interacting with those programs currently in use in the 
CONTROLS package. 


A collection of FORTRAN subroutines, "Optimal Regulator 
Algorithms for the Control of Linear Systems" (ORACLS), by 
E. S. Armstrong of the NASA Langley Research Center, was 
used in the development of the ORACLSX program. Subroutines 
from the "Optimal Systems Control" FORTRAN program (OPTSYS), 
were used as the foundation for the discrete transfer func- 
tion analysis program, TRANFUNC. 

The ORACLS option capabilities include: analog-to- 
digital matrix conversion, transfer function analysis, 
complete eigenvalue analysis, modal analysis, transient 
response calculation, Kalman-Bucy filter synthesis, and 
optimal regulator synthesis. Graphical results are available 
for the transient analysis and transfer function analysis 
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I. INTRODUCTION 


The development of "ORACLS", an interactive capability 
for the study, design , and analysis of digitally controlled 
systems on the IBM 3033, is the subject of this thesis. It 
gives "CONTROLS", a modern control analysis package, the 
ability to handle discrete systems. To that end the FORTRAN 
programs ORACLSX and TRANFUNC were designed. Their develop- 


ment and capabilities are the main focus of this thesis. 


A. BACKGROUND 


The completion of OPTSYSX by Hoden [Ref. 3], made an 
interactive multiple input/multiple output design and anal- 
ysis program available for handling continuous systems. 
OPTSYSX was extended to graphics capable by the programs 
OPTGRAPH and OPTPLOT, by Laptas [Ref. 4] and Diel [Ref. 5], 
respectively. These programs when combined with OPTCALC by 
Diel [Ref. 5] and DACSAP by Cooksey [Ref. 6], became the 
CONTROLS analysis package. They run interactively in the 
VM/CMS environment of the IBM 3033. CONTROLS thus 
constructed, lacked ability to handle discrete systems. 

Optimal Regulator Algorithms for the Control of Linear 
Systems (ORACLS) [Ref. 1] is a collection of FORTRAN subrou- 
tines by E.S. Armstrong of the NASA Langley Research Center. 
It was developed for, "constructing controllers and optimal 
filters for systems modeled by linear time-invariant differ- 
ential or differance equations" [Ref. 1] and employes the 
linear-quadratic-Gaussian (LQG) techniques of modern control 
analysis. The ORACLS package has the capability to handle 
discrete system analysis. A drawback with this analysis 


package however, is the need for the user to code a "driver" 
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program for each problem and to format a data file for each 


run. 


B. PROBLEM OBJECTIVE 


To overcome these problems and fill  CONTROLS' discrete 
analysis void, a program was desired which would use ORACLS' 
discrete system capabilities and interact with those 
Programs already in CONTROLS. 

ORACLSX and TRANFUNC were designed and coded to meet 
those requirements. They can handle various design and 
analysis problems germane to modern control theory. The 
continuous system output from OPTSYSX can be used as input 
to ORACLSX for conversion to a discrete system. The output 
of ORACLSX can be used as input to either of the graphing 
programs OPTGRAPH or OPTPLOT. It can also be used as input 
to the discrete analysis program TRANFUNC. 

It is assumed that the reader is familiar with sampling 
theory and the basic concepts of optimal control theory as 
applied to discrete systems. Wherever possible in this 
thesis and the FORTRAN programs, the nomenclature of ORACLS 
by Armstrong [Ref. 1] was maintained. A glossary of abbrevi- 


ations and matrix notation is provided as page 9. 


C. PREVIEW 


An overview of the CONTROLS package and some modifica- 
tions to it is presented first. The  "ORACLS" option of 
CONTROLS will then be discussed. This encompasses the 
ORACLS EXEC with its program options апа other available 
features. Discussion of the ORACLSX FORTRAN program follows 
next. Each option will be covered individually and its capa- 
bilities demonstrated with an example problem. The final 
discussion will cover the TRANFUNC FORTRAN program. An 


example problem will be used to show its capabilities. 


P2 


II. THE CONTROLS ANALYSIS PACKAGE 


A. GENERAL BACKGROUND 


"CONTROLS" is a modern control analysis package avail- 
able to the system user on the ІВМ 3033 at the Naval 
Postgraduate School, Monterey, California. It is accessible 
in the VM/CMS environment and requires only that the user's 
available storage be increased to the 1 mega byte (1M) 
limit. The system has a 2M capability for running very large 
systems where the capability of the 1M storage limit is 
exceeded. The analysis package may be accessed by simply 
typing; "CONTROLS" on any of the IBM 3278 terminals after 
having defined storage as IM. 


B. THE MANAGE EXEC (CONTROLS' OPTIONS) 


The Manage Exec is used to control the program options 
and features available in the CONTROLS package. Figure 2.1 


is a copy of the controlling menu. 
li) HELP 


The Manage Exec's Help option is a brief description 
of program names and capabilities. See figure 2.2 for the 
HELP menu. 


2. Example Problem 


Option number 5 of the Manage EXEC allows the user 
to obtain a sample OPTMAT data file and a documentation file 
to enter the OPTSYS or ORACLS program options. The instruc- 
tion files OPTSAM and ORCSAM explain the program options of 
the respective OPTSYSX and ORACLSX FORTRAN programs. They 


provide the user a step by step guide to running the 


ЕЗ 
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Figure 2.1 CONTROLS Option Menu. 


programs with the example data and also provide program 


results for comparison. 
3. Program Documentation 


The Program Documentation option lists for the user 
the Master's Theses which generated the programs available 
in CONTROLS. These are the theses as listed in references 
go 


C. MODIFICATIONS TO THE CONTROLS PACKAGE 


This section details modifications made to the OPTSYS 
EXEC and the FORTRAN programs of the CONTROLS package. Some 
of these modifications were necessary to provide the proper 
interface for the addition of the ORACLS option. Some were 
programing carryovers from the ORACLSX program. All FORTRAN 
programs have been updated to VS FORTRAN Release 4.0. This 
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DACSAP: 
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DISCRETE SYSTEM ANALYSIS: 
ANALOG TO DIGITAL MATRIX CONVERSION, KALMAN-BUCY 
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Figure 2.2 CONTROLS HELP Menu. 


update allows the use of designated "character" variables 
and the Block IF statement with the IF THEN, ELSE IF, ELSE 
structure. These changes facilitate programing allowing the 


removal of many confusing "GO TO" statements. 
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1. OPTSYS EXEC Modifications 
a. Data File Protection Checks 


To ensure that the user could not input a 
discrete system (OPTMATD) data file in the OPTSYSX program, 
a flag (IDOPTD) signaling the type of data, was added to the 
OPTMAT and OPTMATD data files. The OPTSYS EXEC was modified 
to check this flag for correctness. Any improper data file 


will be rejected and an error message issued. 
b. Output Selection Option 


Occasionally a user may like to make multiple 
runs through the OPTSYSX program saving the results from 
each run in a different file. Àn option was added to allow 
the selection of six different output files. The user 
selects to either view the program results on the screen or 


to have them sent to the designated listing file. 
c. SHERPA Printer 


The graphing programs OPTPLOT and OPTGRAPH 
generate a DISSPLA METAFILE which through DISSPOP can be 
sent to the SHERPA laser printer. An option was added to the 
OPTSYS EXEC which allows the user to generate thesis accep- 
table plots on the SHERPA printer. 


2. ОРТ5Ү5Х Program Modifications 


Program modifications to some of the subroutines 
Supporting the OPTSYSX program were necessary to ensure a 
proper interface with the ORACLSX program. ORACLSX was 
developed using VS FORTRAN and some of the programing ideas 


from it were carried over to the OPTSYSX program. 
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а. OPTMAT Data File 


AS previously mentioned the IDOPTD flag was 
developed to signal whether the data was continuous or 
discrete. IDOPTD is written as "0" to the OPTMAT data file 
by the subroutine WRTMAT. It is written as a "l" to the 
OPTMATD data file by the subroutine WRMATD. The value of 
IDOPTD is then checked by the OPTSYS AND ORACLS  EXECs, 
allowing only input of proper data type. An error message is 
generated if a conflict arises. 

The OPTMAT data file contains the system parame- 
ters and program matrices. Two changes were made to the 
subroutines RDMAT and WRTMAT to read and write, respectively 
from it. Originally all the matrices were written to the 
file as a single group of numbers making it hard to deci- 
pher. The first change added the corresponding matrix names 
to the file, separating the matrices for ease in identifica- 
tion. The second modification changed the two subroutines, 
writing and reading to and from the file only those matrices 
with non-zero dimensions. This modification prevents the 
Program from writing undefined numbers to the data file 
which are read later as program underflows and cause error 


messages. 
b. ORACLSX Carryovers 


OPTSYSX originally had 14 different subroutines 
to read in the matrices and vectors. These 14 subroutine 
were replaced by a single subroutine READMX. READMX was 
modeled on the subroutine READMD from ORACLSX. It uses the 
matrix identification flag IMAT and the IF THEN, ELSE IF, 
ELSE structure of VS FORTRAN to select the proper name and 


matrix dimensions for the matrix desired. 
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3. OPTCALC and OPTPLOT Program Modifications 


The program results from ORACLSX are in discrete 
time and those from OPTSYSX in continuous time. The OPTPLOT 
program has the capability to plot the results from both 
programs. It was desired to maintain the discreteness of the 
ORACLSX data even when plotting. A data flag, IPLOT, was 
written to the OPTPLOT data file by the OPTCALC and ORACLSX 
programs. A "0" for IPLOT indicates the data file is from 
OPTCALC and should be graphed in the continuous sense. А "1" 
for IPLOT indicates the data file is from the ORTRAL portion 
of ORACLSX and should be graphed in the discrete sense. The 
OPTPLOT program reads the IPLOT flag from the OPTPLOT data 
file. It plots the continuous results connecting the data 
points with а smooth line. For the discrete results it 
plots individual unconnected markers. The markers are spaced 


at integer intervals of the sampling time. 
4.  OPTGRAPH Program Modifications 


TRANFUNC does transfer function analysis of discrete 
Systems in the Z domain.  OPTSYSX does its transfer function 
analysis in the S domain. The OPTGRAPH program can graph the 
results from either program. Data input files for OPTGRAPH 
are: OPTGROL, ОРТСЕМО, апа ОРТСЕСМ. А modification was 
required to indicate their source program. Use was made of a 
data value already available in the data files to indicate 
the source program. A "O" indicating the OPTSYSX program 
and a "1" the TRANFUNC program. The data value is read by 
the OPTGRAPH program and passed to the subroutine RTLO to 


indicate whether a Z-plane template is required on the root- 


locus plot. The template is a plot of lines of constant 
natural frequency and constant damping. See Figure 5.4 The 
template covers the upper half of the unit circle. Its zero 


damping ratio line serving as the boundary between stable 
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III. THE ORACLS EXEC (ORACLS' OPTIONS) 


А. BACKGROUND 


The ORACLS EXEC is used to control the ORACLS option 
available in CONTROLS. It is written іп the CMS EXEC and 
EXEC 2 languages [Ref. 8]. It was modeled after the OPTSYS 
EXEC and includes features similar to those available there. 
Flow in the ORACLS EXEC parallels that in the OPTSYS EXEC, 
allowing the OPTSYS user immediate familiarity with its 
operation. The EXEC sets up the OPTMAT and OPTMATD data 
files. It also performs the data file check as will be 
explained in Section C to follow. The selection of a FORTRAN 
Program option causes the EXEC to set the proper filedefs 
and call the program module. Program execution follows 


Shortly thereafter. ! 


B. AVAILABLE OPTIONS 


1. Program Options 


A general explanation of the capabilities of the 
program options is included in Figure 3.1 An explanation of 
the program and data file connections can be found in Figure 
3.4 of the HELP section to follow. 


2. Output Selection 


Figure 3.2 is a copy of the output selection menu 


from the  ORACLS EXEC. It allows the user the option of 
sending the program results of the ORACLSK and TRANFUNC 
programs to 6 different listing files. Upon exiting 


CONTROLS, the listing file containing the input matrices and 
results, may be printed. The user may select to have the 


output redirected to the terminal Screen at any time. 
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**** THE ORACLS EXEC ALLOWS THE FOLLOWING OPTIONS: 


1 ORACLSX FORTRAN: DISCRETE SYSTEM ANALYSIS 
A-D CONVERSION, OPT. KALMAN- 
BUCY FILTER TAND TRANSIENT ANALYSIS) 


TRANFUNC FORTRAN: -DISCRETE TRANSFER FUNCTION ANAL.- 
COMPENSATOR, 
TRANSFER 
FUNCTIONS, MODAL MATRICES ) 


3. OPTPLOT FORTRAN: (PLOTTING OF TRANSIENT RESPONSE'S 
STATES AND CONTROLS) 


4 OPTGRAPH FORTRAN: (PLOTS: POLE-ZERO, ROOT-LOCUS) 

э ОЮ ТРЫТ SELECTION -- TERMINAL OR "A" DISK FILE 

6 HELP: PROGRAM & DATA FILE RELATIONSHIP EXPLANATIONS 
ТОО ЕК RETURN TO CMS 

-SELECT AN OPTION, ANY OTHER INPUT RETURNS YOU TO MENU- 





Figure 3.1 ORACLS EXEC Options. 


exe OUTPUT OPTION "jx 
"ORACLS LISTING" OR "TRANFUN LISTING" FILES, 
CONTAINING THE P ELE PROGRAM RESULTS OF 6 
DIFFERENT RUNS MAY BE WRITTEN TO YOUR A’ DISK. 
THE OUTPUT FILE NAME IS SELECTED FROM THE FOLLOWING: 


ORACLSI LISTING 


LISTING 
LISTING 
LISTING 
LISTING 


N3 ISTING Al 
DIRECT OUTPUT TO TERMINAL 


ENTER: 1, 2, 3 5 | а, FOR YOUR, OUTPUT FILE NAME 
R ' FOR TERMINATD OUTPU 





Figure 3.2 ORACLS EXEC Output Options. 
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JT HEEE 


The ORACLS EXEC contains a "HELP" option. This 
option explains the data files and their program connec- 
tions. It consists of three separate screens. HELP's screen 
"A" (Figure 3.3) shows the programs available giving a quick 
synopsis of each. It also lists the input and output files 
associated with each program. Screen "B" (Figure 3.4) is a 
list of the data files used in OPTSYS and ORACLS. Screen "C" 
(Figure 3.5) is used to inform the user of disk space limi- 


tations that may arise when using the plotting programs. 


A) PROGRAMS: 


"ORACLSX" FORTRAN IS THE PRIMARY ае РКОСКАМ. 
USES OPTMAT AND OPTMATD. GENERATES OPTMATD. 
ITS ORTRAL PORTION GENERATES THE OPTPLOT 
DATA FILE FOR THE OPTPLOT PROGRAM. 


"TRANFUNC'" FORTRAN DOES TRANSFER FUNCTION ANALYSIS. 
IT GENERATES THE OPTGRAPH FILES: OPTGROL, 
OPTGRNO, AND OPTGRCM FOR OPTGRAPH PROGRAM. 


"OPTPLOT" FORTRAN PLOTS THE TRANSIENT ANALYSIS 
DEVELOPED IN THE ORTRAL PORTION OF ORACLS. 
USES OPTPLOT DATA FILLE FOR PLOTTING ON ТЕК 
618 OR VERSATEC. THESIS ACCEPTABLE PLOTS 
AVAILABLE VIA SHERPA PRINTER. 


"ОРТСВАРН"  FORTRAN PLOTS TRANSFER FUNCTION ANALYSIS 
DEVELOPED IN TRANFUNC PORTION OF ORACLS, 
USES "OPTGROL", "ОРТСЕМО", AND "OPTGRCM 
БАТА FILES FROM TRANFUNC FOR PLOTTING ON 
TEK618, VERSATEC PLOTTER, SHERPA PRINTER. 





Figure 3.3 HELP Screen "A". 
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В) DATA FILES: 


OPTMAT: IS GENERATED BY OPTSYSX AND CONTAINS ALL INPUT 
MATRICES AND REGULATOR AND ESTIMATOR FEEDBACK 
GAINS. CÀN BE USED AS AN INPUT TO THE ORCONV 
SECTION OF ORACLS. 


OPTMATD: IS GENERATED BY ORACLS AND CONTAINS ALL INPUT 
MATRICES AND REGULATOR AND ESTIMATOR FEEDBACK 
mee ID CAN BE USED TO RERUN PROBLEMS 


OPTPLOT: GENERATED BY TRANFUNC AND CONTAINS THE EIGEN- 


VALUE AND TRANSFER FUNCTION ANALYSIS RESULTS. 
T. CAINS, P ZEROS, RESIDUES AND EIGEN 

VALUES IT IS USED BY OPTPLOT FORTRAN 

PLOTTING ON THE TEK618, VERSATEC, OR SHERPA. 


OPTGROL: OPEN LOOP TRANSFER FUNCTION 
OPTGRNO: NOISE TRANSFER FUNCTIONS THRU CLOSED LOOP SYS. 
OPTGRCM: COMPENSATOR TRANSFER FUNCTION 


THE ABOVE THREE FILES ARE GENERATED BY TRANFUNC 
AND USED BY OPTGRAPH TO DO CLASSICAL ANALYSIS. 
PLOTTING VIA THE TEK618, VERSATEC, OR SHERPA 





Figure 3.4 HELP Screen "В". 


C) NOTES: 


IF THE REQUIRED DATA SET IS ON YOUR "A" DISK A 
PROGRAM MAY BE RUN DIRECTLY. 


ALL PEEING PROGRAMS USE DISSPLA AND MAY RE E 


UP TO 2075 OF YOUR "A" DISK SPACE. MAKE SURE 
SPACE IS AVAILABLE FOR THE DISSPLA METAFILES OR 
ЕЕ МЕ55АСЕ STATING THAT YOUR DISK IS FULL 


У ANY INPUT RETURNS TO MENU 





Figure 3.5 HELP Screen "С". 
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C. THE DATA FILES 


Various data files аге generated by the programs 
contained in CONTROLS. This section looks at the the use 


and control of those files. 


1.  OPTMAT and OPTMATD Data Files 


— A A они n  —-— P 


ОРТМАТ апа ОРТМАТП аге data control files. They 
contain the matrices as input by the user or as generated by 
the program. The OPTMAT and OPTMATD files allow the user to 
rerun a problem without having to reenter the matrices each 
time. These files are originally set up on the user's "A" 
disk by the ORACLS EXEC. If selected at the end of each 
program run, the ORACLSX program will write the system 


parameters and matrices to the OPTMATD file. In the begin- 
ning of the next program run, the user will be asked if he 
wants to use some, all, or none of the matrices that were 
saved. 


2. Program Links 


The OPTMAT, OPTMATD, and other data files serve as 
the links between the various programs in CONTROLS. The 
OPTMAT data file is the connection between the OPTSYSX and 
ORACLSX programs. It can be used as input for the ORCONV, 
analog-digital conversion option of ORACLSX. The OPTMATD 
data file is the link between all of the options in ORACLSX. 
The output from ORCONV is written to the OPTMATD file and 
can then be used as the input for any of the other available 
options. The OPTPLOT and OPTGRAPH data files are used to 
link the ORACLSX program to the graphing programs OPTPLOT 
AND OPTGRAPH. 
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3. Protection Features 


The data file protection checks, discussed in the 
modifications to the OPTSYS ЕХЕС, were also incorporated 
into the ORACLS EXEC. These checks prevent an improperly 


renamed data file from being used. See Chapter II, 
С. 


Section 
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IV. THE ORACLSX FORTRAN PROGRAM 


A. PROGRAM OVERVIEW 


ORACLSX is an interactive extension of the discrete 
system capabilities of Optimal Regulator Algorithms for the 
Control of Linear Systems  (ORACLS) [Ref. 1]. It features 
modular programing and double precision computations. Its 
main program iS menu driven and controls four major analysis 
subroutines. In addition to the ORACLS subroutines it has 
numerous subroutines for handling matrices and system param- 
eters. The matrix handling subroutines are modifications of 
those used in OPTSYSX. 

ORACLSX operates under the VS FORTRAN programing 
language. It makes extensive use of VS's character variables 
and Block IF statements. The IF THEN, ELSE IF, ELSE struc- 
ture of the Block IF statement facilitates programing, 
allowing fewer confusing "GOTO" statements, and in fact 
fewer subroutines. See Section D.2. 

The program has been dimensioned to handle systems with 
up to 32 states, 10 controls, 10 measurements/observations, 
and 10 noise inputs; (32,10,10,10). Larger systems, up to 
(0873053053009 can be handled by invoking the alternate 
dimension statements provided in the main program and major 
analysis subroutines. Use of the large system capabilities 
may exceed the storage limits of the ІМ system. If so, the 


2M machine must be accessed. 


B. INCORPORATION OF ORACLS 


The ORACLS package subroutines were set up on the system 
and updated to VS FORTRAN. The update involved removing the 


use of the Hollerith field from all statements that were not 
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МЕ спе Hollerith field from all statementsmthat were not 
format statements. This change affected all the program 
calls to the subroutine PRNT [Ref. 1: pp. 8,9]. REAL decla- 
rations were replaced by CHARACTER declarations for appro- 
priate variables in the subroutines RDTITL [Ref. 1: p. 5], 
and PRNT. The "IMPLICIT REAL*8  (A-H,0-Z)" statement was 
added to all subroutines, updating the package to double 


precision. Integer variables declared as REAL, were redec- 
lared as  REAL*8. Double precision function names were 
substituted for "ABS",  "ALOG10", and "AMAX1" [Ref. 9: pp. 


303,304]. After completion of the conversion to VS FORTRAN, 
the four example problems included in the ORACLS manual 
[Ref. 1], were run. Program results generated by the updated 


package agreed completely with those given in the manual. 


С.  ORACLSX'S MATRIX NOMENCLATURE 


The discrete system nomenclature of ORACLS [Ref. 1] was 
maintained wherever possible in ORACLSX. Exceptions are the 
(GAMD), (V1), and (V2) matrices. {A} is the System, Plant, 
or State Weighting Matrix. {B} is the Control Weighting 
Matrix. САМО) is the Process Noise Weighting Matrix. The 
{H} Matrix is the Weighting of the Observations. {G} is the 
Measurement-Feed Forward Distribution Matrix. (F) is the 
Control Gain Matrix. (K) is the Filter Gain Matrix. (Q) is 
the Output Cost Matrix. (R) is the Control Cost Hatrix. (V1) 
is the Process Noise Intensity (PSD) Matrix. (V2) is the 
Measurement Noise Intensity (PSD) Matrix. and (S) is the 


Cross Product Weighting Matrix. 


D. SYSTEM/MODEL DESCRIPTION 


ORACLSX works with discrete time systems. It bases its 
analysis on difference equations. In the following differ- 


ence equations, "T" is the Sample Time Interval. 
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The SYSTEM Equation: 


X((I*l)*T) » (AJX(I*T) * (BJU(I*T) * (GAMDOWD(I*T) (4.1) 


The OUTPUT Equation: 


Ү(Т) = (HIK O CUT) (2 


The MEASUREMENT Equation: 


Z(I) Ey ТЕП (a) (4.3) 


The CONTROL LAW: 


U( 1) = state PX) (4.4) 
Where: 
X(I) - State Vector 
Y(I) - Output Vector 
Z(I) - Measurement Vector 
WD(I) - White Process Noise Vector 
V(I) - White Measurement Noise Vector 


The "WD" vector has zero mean value, and a covariance matrix 
{Vl}, where 


E(WD) = O (4.5) 
and 


(Vl) » E(WD*WD!) (4.6) 
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The "V" vector has zero mean value, and a covariance matrix 
(V2), where 


E(V)= O (4.7) 
апа 
{V2} = E(VeV") (4.8) 
The quadratic performance (or Cost) index for the linear 


quadratic regulator is the expected value of: 
The COST EQUATION: 


LEM lim Z ( X" (I*T)(VI)EK(I*T) 9 X" (I*T)C)X(I*T)... 
..+ UT(IsT)(V2)U(I1ST) ) (4.9) 


E. MAIN PROGRAM 


1. Error Suppression 


VS FORTRAN has an extensive error message library. 
Program execution is terminated when the number of error 
messages generated exceeds the default value. In ORACLSX 
calls to ERRSET are made to allow unlimited "Overflow", 
"Underflow'", and "Divide Check" error messages. When these 
errors occur the standard corrective fix-up takes place 
allowing program continuation. ERRSET calls were also made 
to prevent "Illegal Decimal Character" and "Dimension Check" 


error messages from terminating the program. 
2. DATA Handling 


Data input into the ORACLSX program can be done 
interactively, one matrix at a time, or all at once, using 
an OPTMAT/OPTMATD file. 
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а. System Parameters and Save Flags 


If an old OPTMAT/OPTMATD file does not exist the 
program will prompt the user for the system parameters. The 
system parameters (number of states, number of controls, 
number of observations/measurements, number of process noise 
inputs) define the size of the system to be analyzed. If 
there is an OPTMAT/OPTMATD file available, the user will be 
given a choice as to how many of the old matrices are to be 
Saved and used again. His choices are: 1) Use all of the old 


matrices, 2) Use some of the old matrices, 3) Input all new 


matrices. If the user selects to use none of the old 
matrices, the program will prompt for all of the system 
parameters. If he selects to use some or all of the old 


matrices, save flags will be set for those matrices that are 
to be saved. The program will now prompt for those system 


parameters that were previously zero and are now required. 
b. OPTMAT AND OPTMATD DATA Files 


ORACLSX can input system parameters and matrices 
from either the OPTMAT or OPTMATD data files. OPTMAT is used 
as input to ORACLSX's  ORCONV option. OPTMATD is used as 
input for all other options. System parameters and matrices 
are read from the OPTMAT and OPTMATD files using the RDMAT 
and RDMATD subroutines, respectively. An OPTMAT file may be 
written by the WRTMAT subroutine if the user wants to save 
the continuous matrices input into the ORCONV option. An 
OPTMATD file may be written by the WRMATD subroutine at the 
end of all program options. The OPTMATD file will contain 
the discrete matrices generated by the program during that 
run as well as those matrices input by the user. 

The OPTMAT file can be used for reruns through 
the ORCONV option or as input to the OPTSYSX program. The 
OPTMATD file can be used as input to any of ORACLSX's other 
options or as input to the TRANFUNC program. 
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ОО АШ ПГегасгтуе Input 


The subroutines READMX and READMD were developed to 
handle the input of the continuous and discrete matrices, 
respectively. Format statements for the desired matrix are 
selected using a matrix identification flag, IMAT, submitted 
to the subroutines. 

The NULL, UNITY, and SCALE subroutines of ORACLS 
were adapted to READMX and READMD to have the program 
generate the corresponding matrix. This option relieves the 
user of the time consuming task of entering certain matrices 
an element at a time. For identity matrices or those with a 
large number of zero elements this proves to be a valuable 
feature. The generation of unity and unity scaled matrices 
1s only available if the matrices are square. 

The matrix input sections of READMX and READMD are 
modeled after those originally used in OPTSYSX. As in 
ОРТ5У5Х the user has the ability to change any matrix 


element before program execution. 


4. Program Required User Inputs 


Three subroutines; RDREAL, RDINT, and RDCHAR were 
borrowed from OPTSYSX [Ref. 3: р. 20] to read in program 
required user inputs. Subroutine RDREAL is call to input the 
user's response at points where a real number or Zero 
integer may be expected. Subroutine RDINT is called to input 
the user's response at any point where a non-zero integer 


may be expected. It is used to input all program option 


selections. Subroutine RDCHAR is used whenever a "(Y)es” or 
"(N)O" answer is required. RDCHAR was updated in ORACLSX to 
use the character variable available in VS FORTRAN. In 
ORACLSK and TRANFUNC, RDCHAR is used to read the first 
letter of the user's input. The user may enter "Y" or 
"YES", or "N" or "NO" as an acceptable responses. 
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In each of these three subroutines a protection 
feature is built in. The user may enter one "null" entry 
without program execution  halting. If this happens, a 
warning message will be sent to the screen and the program 
will recover. Two "null" entrees in a row will produce an 
error message and terminate the program sending the user 
back to the controlling EXEC. This is a handy feature as it 
allows an escape from the program anywhere the user is 
prompted for an input. Also in these subroutines is a check 
for correct data type. If the data input by the user is not 
of the correct type, an error message is sent to the screen 
instructing the user to enter the correct type (i.e. real, 


integer, character). 


F. PROGRAM CAPABILITIES 


This section explains the analysis capabilities of the 
ORACLSX program. Figures 4.1 and 4.2 show ORACLS's main and 
Discrete System Analysis (ORDSAL) option menus. 


ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
ORDSAL: DISCRETE SYSTEM ANALYSIS 

ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
HELP: PROGRAM & EQUATION DESCRIPTIONS 


EXIT: RETURN TO ORACLS EXEC 


SELECT OPTION - 5; ТЕ YOU ARE USING THE OUTPUT 
TO DISK OPTION AND WANT TO CHANGE FILENAMES. 
NORMALLY ANY TIME DURING THE PROGRAM TWO 
SUCCESSIVE NULL "ENTERS" TERMINATES THE PROGRAM 
SENDING YOU BACK TO THE ORACLS EXEC. 


SELECT AN OPTION: 1,2, 3, 4, OR 5. 





Figure 4.1 ORACLSX's Main Option Menu. 
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ORDSAL OPTIONS: 
KBFIL : DISCRETE KALMAN-BUCY FILTER ANALYSIS 
OPTREG: DISCRETE OPTIMAL REGULATOR ANALYSIS 
EXIT: RETURN TO ORACLS OPTIONS 


EXIT: RETURN TO ORACLS EXEC 
SELECT AN OPTION: 1, 2, 3, OR 4. 
NORMALLY ANY TIME DURING ANY PROGRAM 


TWO 
SUCCESSIVE NULL "ENTERS" TERMINATES THE PROGRAM 
SENDING YOU BACK TO THE ORACLS EXEC. 





Figure 4.2 Discrete System Analysis Option Menu. 


For each of the program options a recorded terminal 
session demonstrating that option is available as a separate 
appendix. The HELP option shown in Figure 4.1 will be 
discussed in Section F of this chapter. 

The discrete analysis subroutines of ORACLS are used as 
the foundations for the options in ORACLSX. In the discus- 
sions to follow references will be made to the ORACLS manual 
[Ref. 1] where an in depth discussion of each subroutine is 


provided and original references cited therein. 
1. ОКСОМУ Option 
а. Program Operation 


The purpose of the ORCONV option of ORACLSX is 
to perform analog to digital matrix conversion. A user 
entering ORACLSX with the results from the OPTSYSX program 
or other continuous matrices must first convert the input 
data to discrete values. With the discrete matrices from 
ORCONV the rest of the ORACLSX options maybe exercised. 


9IS 


ORCONV conversion capabilities are as shown іп 
Figure 4.3 Those listed in the OPTSYS column are continuous 
matrices and maybe converted to the discrete matrices 
listed in the ORACLS column. A discrete cross product 
weighting matrix, {S}, will be generated if the {A} and {B} 
continuous matrices are converted to the discrete {Q} and 


{R} matrices. 


MATRIX DESCRIPTION OPTSYS NAME ORACLS NAME 


SYSTEM MATRIX: 
CONTROL MATRIX: 


NOISE MATRIX: 

OUTPUT COST MATRIX: 

CONTROL COST MATRIX: 

PROCESS NOLSE MATRIX: 
MEASUREMENT NOISE MATRIX: 

CROSS PRODUCT WEIGHTING MATRIX: 





Figure 4.3 ORCONV Conversion Capabilities. 


b. Input and Output Data Files 


ORCONV will accept input from the screen or 
allow the use of an old OPTMAT data file. The OPTMAT file 
may be from a previous ORCONV run or one from the OPTSYSX 
program.  ORCONV outputs the discrete matrices generated to 
an OPTMATD data file which may be used in the other ORACLSX 
options or the TRANFUNC program. 


с. Solution Algorithms 


Conversion of the matrices is carried out using 
the subroutines EXPINT [Ref. 1: pp. 49-50] and SAMPL 
[Ref. 1: pp. 61-63]. EXPINT is called for conversion of the 
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(F) matrix. Its results are also used for conversion of the 
(G) matrix. SAMPL is called to compute the (Q), (R), and (S) 
matrices. ТЕ is also called to compute the {Vl} and (V2) 
matrices. 

EXPINT's computations are done using finite- 
series algorithms [Ref. 1: p. 49]. SAMPL uses another set 
difinite series algorithms [Ref. 1: p. 61] for its computa- 
tions. Both subroutines can be effectively thought of as 
part of a sample and hold device. Figure 4.4 shows ап 
example of how the A/D conversion may be done using a zero- 


order hold device and a sampler. 






Zero-order hold Sampler 


Figure 4.4 Model of the A/D Converter. 


A signal is fed to the zero-order hold, (ZOH) 
device which "holds" the input signal for atime period 
equal to the sample time interval. From the ZOH the signal 
passes through the analog system to the sampler. Here the 


sampler "samples" the signal at the prescribed time interval 


"T" to produce a discrete signal. For ORCONV the sample 
time interval "T", is input by the user. Some texts refer to 
the sample time interval as "delta" or "delt". Іп this 
program in will be called,  "delt" and in equations listed 


will be represented by "T". 
а. ORCONV HELP 


Та the HELP option of  ORAGLSX, a section 
covering the ORCONV option is available. Information avail- 


able there is included in Appendix A; "ORCONV Example Run". 
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e. Example Problem Run 


To demonstrate the capabilities of the ORACLSX 
program an example Problem from Linear Optimal Control 
Systems by Kwakernaak and Sivan [Ref. 2]. was used. The 
continuous-time positioning system of Example 2.4 [Ref. 2: 
pp. 133-136] is updated to a digitally controlled system in 
Example 6.2; ‘Digital Positioning  System' [Ref. 2: РР. 
447-4491. The continuous-time matrices are converted to 
their discrete-time equivalents by the  ORCONV option of 
ORACLSX. 

Appendix A is a recording of the terminal 
session where the example problem data was used to exercise 
the ORCONV option of ORACLSX. The session begins with entry 
in to the CONTROLS analysis package. Here the ORACLS option 


is selected, placing the user in the ORACLS EXEC. From the 
ORACLS EXEC's options the ORACLSX FORTRAN program is 
selected. Before executing the ORCONV option, information 


comparing the OPTSYSX and ORACLSX programs is viewed. HELP 
information on the ORCONV program is also viewed at this 
time including a list of the equations used іп the conver- 
Sion process. 

The terminal session continues as the system 
parameters are input defining the system. Selection of the 
sampling time interval is made anda status screen showing: 
the matrices to be converted, the sample time selected, and 
the system parameters, is generated. Input of the contin- 
uous matrices follows, with the user given an option to save 
them in ап OPTMAT file. The continuous matrices input and 
the discrete matrices generated appear next. These are the 
discrete matrices as shown in equation 6-26 of [Ref. 1: p. 
448]. At this point the user is given the option of saving 
the discrete matrices generated in an OPTMATD file. The 
ORCONV option is now complete and the user is sent back to 


ORACLSX's options from where the program is exited. 
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Za KBFIL Option 
a. Program Operation 


The purpose of KBFIL is to solve the discrete time-invariant 
asymptotic optimal Kalman-Bucy filter problem. The optimal 
filter problem is to construct an estimate of the state, 
X(I*T) from knowledge of the measurements,  Y(I*T) and the 
control inputs, U(I*T) up to the time (I-1)*T. 

Given the estimate of X, denoted as R. the esti- 


mator equation is then written as: 
$(r41) - (A)X(I)+(B)U(I)+(K)(Y(1)-(H)XXI1)-(G)U(I)) (4.10) 


Solution of this problem produces the optimal filter gain 
matrix; (K). Modification of format statements in ASYMFI 
was required. In all output statements referring to the 
filter gain matrix the "F" matrix was renamed the "К" 
matrix. All other output statements were updated to use the 
"V1", "V2", and "GAMD'" nomenclature of ORACLSX. 


b. Input and Output Data Files 


KBFIL will accept input from the screen or allow 
use of an old OPTMATD data file. The OPTMATD file may be 
from a previous KBFIL run or from an  ORCONV run. KBFIL 
outputs the computed optimal filter gain, "К" matrix, апа 
input matrices to an OPTMATD file. This file may then be 
used as input to the OPTREG option which will generate the 
control gain matrix, "Е". With these two gain matrices the 
Compensator analysis option of the TRANFUNC program may be 


exercised. 
c. Solution Algorithms 


Solution of the Kalman-Bucy filter problem is 
carried out using the ORACLS's ASYMFI subroutine. [Ref. 1: 


37 


рр. 94-99|. Duality of the Optimal Observer and the Optimal 
Regulator [Ref. 2: Sec 4.4] is employed in the solution, 
allowing ASYMFI to invoke the subroutine ASYMRE [Ref. 1: pp. 
88-93] to solve the linear optimal regulator problem. In its 
solution of the regulator problem,  ASYMRE calls the subrou- 
tine RICTNW which solves the discrete steady-state Riccati 
equation by the Newton Algorithm [Ref. 1: pp. 84-87]. 


а. KBFIL HELP 


In the HELP option of ORACLSX, a section 
covering the KBFIL option is available. Information avail- 


able there is included in Appendix В; "KBFIL Example Run". 
e. Example Problem Run 


The Kwakernaak and Sivan discrete control 
example problem data, previously discussed in the ORCONV 
option, was used to demonstrate the KBFIL option. Appendix 
Bis a recording of a terminal session demonstrating the 
KBFIL option. It uses input from the OPTMATD data file in 
conjunction with new data for those matrices not previously 
generated. The session starts with the KBFIL HELP informa- 
tion. Menu screens for the Manage EXEC, and ORACLS EXEC are 
not shown in this session as they are duplicates of those in 
the beginning of the ORCONV Example Run, and may be viewed 
in Appendix A. Menu screens’ for exiting the program have 
also been deleted here as they are identical to those in 


Appendix A. 
3. OPTREG Option 
a. Program Operation 


The purpose of OPTREG is to solve the discrete 
time-invariant asymptotic linear optimal regulator problem 


with noise free measurements. Solution to the problem yields 
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the Control Gain Matrix; {F}. If the system is controlable 


and observable, a solution exists and is given by. 


(F} = {({R}+{B}" (P}{B}}7?{B}* (P}{A} (4.11) 
where 
(Р) = (M}*{(P}{M}+{F}* (R}(F}+{H}"(Q}{H)} (4.12) 
with 
(М) = (A)-CB)(CF) (2.19) 
During Program operation the user is asked 
whether the system is stable or unstable. If ап unstable 
response 15 entered, the matrix {A-BF} is evaluated and 
tested for stability relative to Alpha = .9 using the 
subroutine TESTST [Ref. 1: pp. 43-44]. If a stabilizing 
gain is required, it iS computed in the subroutine DSTAB 


[Ref. 1: pp. 70-73]. 

OPTREG is set up to handle a system with ог 
without a discrete cross product weighting matrix; {S}. If a 
cross product matrix 15 present the system will calculate a 
new gain matrix; {FE}, which will eliminate the cross 
product term in the quadratic scalar function. With the gain 
matrix (FE), new "HAT" matrices are computed for use in the 
System equation and quadratic scalar function. These modi- 
fied equations are then used and the control gain matrix; 
(F), is calculated. Also output is the "FHAT" matrix which 
contains the control gains for the "HAT" system. If no cross 
product term exists, then {FHAT} = {F}. 
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b. Input and Output Data Files 


OPTREG will accept input from the screen or 
allow use of an old OPTMATD data file. The OPTMATD file may 
be from a previous KBFIL, OPTREG, or ORCONV run. OPTREG 
outputs the computed optimal filter gain; "F" matrix, апа 
input matrices to an OPTMATD file. This file may then be 
used as input to the transfer function analysis options 
available in the TRANFUNC program. 


c. Solution Algorithms 


The linear optimal regulator problem is solved 
using ORACLS's ASYMRE subroutine ([Ref. 1: рр. 88-931. 
ASYMRE inturn solves the discrete steady-state Riccati equa- 
tion using the Newton algorithm of subroutine RICTNW 
[Ref. 1: pp. 84-87]. 


d. ОРТКЕС НЕГР 


In the HELP option of  ORACLSX, а section 
covering the OPTREG option is available. Information avail- 


able there is included in Appendix C; "OPTREG Example Run". 
e. Example Problem Run 


Kwakernaak and Sivan example 6.14 [Ref. 2: p. 
494] is a continuation of Example 6.2 [Ref. 2], used in the 
ORCONV Example Run, and is used here to demonstrate the 
OPTREG option of ORACLSX. Appendix C is a recording of a 
terminal session demonstrating the OPTREG option. It uses 
input from the OPTMATD data file in conjunction with new 
data for those matrices not previously generated. The 
session starts with the OPTREG HELP information. Menu 
screens for the Manage EXEC, and ORACLS EXEC are not shown 
in this session as they are duplicates of those in the 


beginning of the ORCONV Example Run, апа may be viewed in 
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Appendix À. Menu screens for exiting the program have also 
been deleted here as they are identical to those in Appendix 
A. 


4.  ORTRAL Option 
a. Program Operation 


The ORTRAL option of ORACLSX computes and prints 


the transient response of the time-invariant discrete system 


out to a desired number of stages. One stage is equal to 
the sampling time interval. Given a set of initial condi- 
tions; Х(0), the transient response is calculated using the 
following: 


System Equation: 


X(I+1) = (A)X(I) + (B)U(I) (4.14) 
Output Equation: 
к=н ту чот) (i 
Control Law: 
ШТІ ЕХ NES BD (4.16) 


Where UC is the Control Input (Driving Function) апа may be 
selected as either a Step or Ramp. Selection of the start 
and stop stages for the control input are selected by the 
user. The user is also given the option of individually 
selecting those control variables which he desires to drive. 
Normally five lines of output are generated for each stage 
computed. An option was added to allow the user to select 
only a desired number of stage results for terminal viewing. 


This option prevents the user from having to spend time 
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needlessly clearing the screen when a large number of stages 
is used. This option does not limit the program output to 
the OPTPLOT data file. A transient response plot for the 
entire number of stages is still possible. 

During program operation the user is asked 
whether the system is stable or unstable. Ifa stable 
response is entered, ORTRAL will compute and output the 


steady state value. 
b. Input and Output Data Files 


ORTRAL will accept input from the screen ог 
allow use of an old OPTMATD data file. The OPTMATD file may 
be from a previous OPTREG or ORTRAL run. ORTRAL outputs to 
the OPTMATD data file the system parameters апа input 
matrices. No results generated by the program are output to 


the OPTMATD file. Transient response results of the states 


and controls are output to the OPTPLOT data file. The 
OPTPLOT FORTRAN program, available in the ORACLS EXEC, uses 
this file. It gives the user the option of obtaining the 


results in graphical format. 
c. Solution Algorithms 


Transient response calculations are carried out 
using ORACLS's TRANSI subroutine [Ref. 1: рр.57-5 
Several modifications were necessary to TRANSI to іпсогро- 
rate the different driving functions and selectivity of 
controls to Бе driven. Output format changes were also 
necessary to conform to the ORACLSX nomenclature. The input 
control; UC was added to the list of outputs, to allow 


viewing of the changes in the driving function. 
d. ORTRAL HELP 


In the HELP option of ORACLSX, a section 
covering the ORTRAL option is available. Information avail- 
able there is included in Appendix D; "ORTRAL Example Run". 
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е. Example Problem Run 


Kwakernaak and Sivan example 6.13 (Ref. 2: p. 
490] was used to demonstrate the ORTRAL option. It isa 
continuation of the digital positioning system problem of 
Example 6.2 (Ref. 2], discussed previously in the ORCONV 
example run. 

Appendix D is a recording of a terminal session 
demonstrating the ORTRAL option. It uses input from the 
OPTMATD file inconjunction with control gains input sepa- 
rately. The session starts with the ORTRAL HELP information. 
Menu Screens for the Manage EXEC, and ORACLS EXEC are not 
shown in this session as they are duplicates of those in the 
beginning of the ORCONV Example Run, апа тау be viewed in 
Appendix A. Menu screens for exiting the program have also 
been deleted here as they are identical to those in Appendix 
А. 

Figures 4.5 and 4.6 show ORTRAL's results for 
the example problem as generated Бу the OPTPLOT FORTRAN 
program. They are included for comparison with those of 
IEPure 6.12 [Ref. i: p. 491]. 


G. HELP 


Figure 4.7 shows the HELP feature available in ORACLSX. 
The menu gives the user the option of obtaining hardcopy of 
the information selected for viewing. HELP is meant to 
provide the user with sufficient insight into program opera- 
tion, nomenclature, and solution methods used. It does not 
provide a detailed discussion of the ORACLS subroutines used 
and their solution techniques. Ап interested user should 
obtain GENCES  [Ref. 11 and refer to the "Solution 
Algorithms" discussions available in the earlier portions of 


this chapter. 
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Figure 4.5 ORTRAL Example Run Plot #l. 
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Figure 4.6 ORTRAL Example Run Plot #2. 
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Information available in HELP has been included іп the 
terminal sessions which demonstrate the program options. See 


Appendices A-D. 


"HELP" PROVIDES, AN CE OF THE PROGRAM OPTIONS 
AVAILABLE IN  "ORACLS" AND SOME OF THE EQUATIONS 
USED IN THE CALCULATIONS. 


-- TERMINAL VIEWING OF SELECTED INFORMATION 
TERMINAL VIEWING PLUS HARDCOPY TO YOUR "A" DISK 


-- EXIT: RETURN TO THE ORACLS OPTIONS 
NOTE: OPTION KA WILL PLACE THE FILE: "ORACLS HELP"' 
ON YOUR "A" DISK CONTAINING THÉ SELECTED 
INFORMATION THAT YOU HAVE VIEWED. 


SELECT OPTION: 1, 2 TORE) 





Figure 4.7 ORACLSX HELP Control Menu. 
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V. THE TRANFUNC FORTRAN PROGRAM 


A. PROGRAM OVERVIEW 


TRANFUNC 15 ап interactive discrete transfer function 
analysis program. It features modular programing and double 
precision computations. Its main program is menu driven and 
controls four major analysis options. Additional options 
provide the user with HELP information or allow the MARKOV 
parameter to be changed (See Section E.5). 

TRANFUNC operates under the VS FORTRAN programing 
language. It makes extensive use of VS's "Character" vari- 
able and Block IF statements. Program dimensions are 
similar to those in ORACLSX and have been set to handle 
systems with up to 32 states, 10 controls, 10 measurements/ 
observations, and 10 noise inputs. For information on its 
capability to handle larger systems, see Chapter IV, Section 
А. 


B. INCORPORATION OF ORACLS AND OPTSYS SUBROUTINES 


System parameter and matrix handling subroutines used in 
TRANFUNC are from ORACLS and ORACLSX. Updating was required 
for the ORACLSX subroutines;  SYSPAR, RMATFD, and INPUT, to 
handle the TRANFUNC options. Subroutines used in the 
transfer function analysis are from OPTSYSX. Formating 
changes were required in subroutines CNORM and MODE to main- 


tain the discrete nomenclature used in ORACLSX. 


С. TRANFUNC'S NOMENCLATURE 


The discrete system nomenclature of ORACLS [Ref. 1] was 


maintained wherever possible in TRANFUNC. An exception was 
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the {GAMD} matrix. (A) is the System; VPE or State 
Weighting Matrix. (B) is the Control Weighting Matrix. 
{GAMD} is the Process Noise Weighting Matrix. The {H} Matrix 
is the Weighting of the Observations. {G} is the 
Measurement-Feed Forward Distribution Matrix. {F} is the 


Control Gain Matrix. {K} is the Filter Gain Matrix. 


D. SYSTEM/TRANSFER FUNCTION DESCRIPTIONS 


TRANFUNC works with discrete time systems. It bases its 
analysis on difference equations. In the following differ- 


ence equations, "T" is the Sample Time Interval. 
The SYSTEM Equation: 
X((I*1)*T) = (A)JX(I*T) * (BJU(I*T) + (GAMD)WD(I*T) (5.1) 
The OUTPUT Equation: 
Y(I) = (H)X(I) * (G)U(I) (57 20) 
The MEASUREMENT Equation: 
Z(I) = (H)X(I) * (G)U(I) * V(I) (5.3) 
Where: 


X(1) - State Vector 
Y(I) - Output Vector 


Z(I) - Measurement Vector 
WD(I) - White Process Noise Vector 
V(I) - White Measurement Noise Vector 


For more discussion of the equations and variables see 
Chapter 4, Section D. 
TRANFUNC does its transfer function analysis in the z 


domain. Its transfer functions are as follows: 
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Open Loop Transfer Function 


ЧИЕН СНО ани ори (5.4) 
Noise Transfer Function 
РСТ ГТ SPS (B (PPS GAMD) бол 
Compensator Transfer Function 
ОИСИ 2 У СА) СВОЕ КУОНУ ^ CK) (5.6) 


E. MAIN PROGRAM 


p Error Suppression 


VS FORTRAN has an extensive error message library. 
Program execution is terminated when the number of error 
messages generated exceeds the default value. In TRANFUNC 
calls to ERRSET are made to allow unlimited "Overflow", 
"Underflow", and "Divide Check" error messages. When these 
errors occur the standard corrective fix-up takes place 
allowing program continuation. ERRSET calls were also made 
to prevent "Illegal Decimal Character" and "Dimension Check" 


error messages from terminating the program. 
2. DATA Handling 


Data input into the TRANFUNC program can be done 
interactively, one matrix at a time, or all at once, using 
an OPTMATD file. 


a. System Parameters and Save Flags 


ТЕ an old OPTMATD file does not exist the 
program will prompt the user for the system parameters. The 


system parameters (number of states, number of controls, 
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number of observations/measurements, number of process noise 
inputs) define the size of the system to be analyzed. If 
there is an OPTMATD file available, the user will be given a 
choice as to how many of the old matrices are to be saved 
and used again. His choices are: 1) Use all of the old 


matrices, 2) Use some of the old matrices, 3) Input all new 


matrices. If the user selects to use none of the old 
matrices, the program will prompt for all of the system 
parameters. If he selects to use some or all of the old 


matrices, save flags will be set for those matrices that are 
to be saved. The program will now prompt for those system 


parameters that were previously zero and are now required. 
b. The OPTMATD DATA File 


TRANFUNC can input the system parameters and 
matrices from an OPTMATD data file.System parameters and 
matrices are read from the OPTMATD file using the RDMATD 
Subroutine. A OPTMATD file may be written by the WRMATD 
Subroutine at the end of all program options. The OPTMATD 
file will contain only the discrete matrices input by the 
user during that run as all results are sent to the screen/ 
listing file and the OPTGRAPH data files; OPTGROL, OPTGRNO, 
OPTGRCM. 

The OPTMATD file can be used as input to any of 
TRANFUNC's options or as input to the ORACLSX program. The 
OPTGRAPH data files can be input to the OPTGRAPH FORTRAN 
program to obtain pole-zero plots and root-locus’ plots for 


the respective transfer functions. 
3. Interactive Input 


The subroutine READMD, developed for ORACLSX, is 
used to handle the input of the matrices. Format statements 
for the desired matrix are selected using a matrix identifi- 
cation ilag, IMIT, submitted to the subroutine by the 


calling program. 
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The NULL, UNITY, and SCALE subroutines of ORACLS 
were adapted to READMD to have the program generate the 
corresponding matrix. This option relieves the user of the 
time consuming task of entering certain matrices an element 
at a time. For identity matrices or those with a large 
number of zero elements this proves to be a valuable 
feature. The generation of unity and unity scaled matrices 
is only available if the matrices are square. The user is 
given the chance to change any matrix element before program 


execution. 


4. Program Required User Inputs 


Three subroutines; RDREAL, RDINT, and RDCHAR were 
borrowed from OPTSYSX [Ref. 3: p. 20] to read in program 
required user inputs. Subroutine RDREAL is call to input the 
user's response at points where a real number or zero 
integer may be expected. Subroutine RDINT is called to input 
the user's response at any point where a non-zero integer 


may be expected. ТЕ is used to input all program option 


selections.  Subroutine RDCHAR is used whenever a "(Y)es" or 
"(N)O" answer is required. RDCHAR was updated in ORACLSX to 
use the character variable available in VS FORTRAN. In 
ORACLSX and TRANFUNC, RDCHAR is used to read the first 
letter of the user's input. The user may enter "Y" or 
"YES", or "N" or "NO" as an acceptable responses. 


In each of these three subroutines a protection 
feature is built in. The user may enter one "null" entry 
without program execution halting. If this happens, a 
warning message will be sent to the screen and the program 
will recover. Two "null" entrees in a row will produce an 
error message апа terminate the program sending the user 
back to the controlling EXEC. This is a handy feature as it 
allows an escape from the program anywhere the user is 


prompted for an input. Also in these subroutines is a check 


EN 


for correct data type. If the data input by the user is not 
of the correct type, an error message is sent to the screen 
instructing the user to enter the correct type (i.e. real, 


integer, character). 


F. PROGRAM CAPABILITIES 


This section explains the analysis capabilities of 
TRANFUNC. Figure 5.1 shows TRANFUNC's main option menu. 
Discussions of each of the options listed in the menu will 
follow. Discussions will cover; the data files output by the 
different options, the solution algorithms used, and the use 
of an example problem to demonstrate each option. All anal- 
ysis options listed in Figure 5.1 allow input from the 
Screen or from an old OPTMATD data file. The OPTMATD file 
may be from a previous TRANFUNC run or from the ORACLSX 
program. All options allow the user to create an OPTMATD 
data file at the end of program analysis. This file will 
save only the discrete matrices entered into the program. No 
matrices generated in the program are written to the OPTMATD 
file. Results to be used in the graphing options of the 
ORACLS EXEC are output to the OPTGRAPH data files. They will 


be discussed in the individual option discussions. 


1. Open-Loop Eigensystem Analysis, {A} Matrix 


a. Program Operation 


Option l performs eigenvalue analysis on the 
system matrix, {A}. {A} is a Square matrix dimensioned, 
(NS,NS). "М5" eigenvalues will be found and printed out. 
Eigenvalues may be real or complex. The open-loop right 
elgenvector matrix; {T}, and open-loop left eigenvector 


matrix; (T)^!, are also printed out. 
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TRANFUNC OPTIONS: 
OPEN-LOOP EIGENSYSTEM ANALYSIS ONLY, (A) MATRIX 
OPEN-LOOP TRANSFER FUNCTION & EIGENSYSTEM 
NOISE TRANSFER FUNCTION & CLOSED-LOOP EIGENSYSTEM 
COMPENSATOR TRANSFER FUNCTION & ESTIMATOR EIGENSYS 


MARKOV PARAMETER CHANGE 
HELP -- PROGRAM DESCRIPTION 
EXIT -- RETURN TO ORACLS EXEC 
SHMECT OPTION: 1, 2, 3, 4, 5, 6, 7; (SETS ITFOL FLAG) 
NOTE: MODAL ANALYSIS AVAILABLE IN OPTIONS: 2, 3, 4 





Figure 5.1 TRANFUNC Option Menu. 


b. Output Data Files 


This option outputs the system parameters апа 
eigenvalues to the OPTGRAPH file; OPTGROL data. 


c. Solution Algorithms 


Analysis of the eigensystem is carried out using 
the following subroutines from OPTSYSX [Ref. 3]: BALANC, 
ORTHES, ORTRAN, EREXIT, and BALBAK. The OPTSYSX subroutine 
CNORM, updated to the discrete nomenclature and retitled 
DCNORM is also used. 


d. Example Problem Run 


To demonstrate the eigensystem analysis option 
of TRANFUNC, the Kwakernaak and Sivan discrete control 
example problem, discussed in ORACLSX, was used again. 
Appendix E is a recording of the terminal session where the 


example problem data was used to demonstrate the eigensystem 
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analysis option. Also included there is a copy of the HELP 
information available in the TRANFUNC program. 


2. Open-Loop Transfer Function Analysis 


Figure 5.2 shows the Open-Loop Transfer Function 
option menu. It contains the open-loop’ transfer function 
equation and available analysis options. The user is 
informed here that to use the OPTGRAPH plotting program 
requires the selection of either option 1 or 2. Also pointed 
out are the nonzero requirements for the system parameters 
“МС Шапа “О! the number of controls and number of 


observations/measurements respectively. 


OPEN-LOOP TRANSFER FUNCTION OPTIONS: 
A/U = (Ну z(I)- (A) )INV (D) 
OPTION l -- POLES, RESIDUES; ANDITZEROSTCOMEUMTED 
OPTION 2 -- ONLY POLES AND ZEROS COMPUTED 


OPTION 3 -- ONLY POLES AND RESIDUES COMPUTED 


OPTION & -- EXIT -- RETURN TO TRANFUNC OPTIONS 
SELECT OPTION: 1, 2, I FORT: (SETS ITFl=1,2,0R,3) 
NOTE: POLES AND ZEROS MUST BE SELECTED IF YOU 


WISH TO USE OPTGRAPH. 
NOTE: REQUIRES NC>0, NO>0 





Figure 5.2 Open-Loop Transfer Function Options. 


a. Program Operation 


The user selects the desired option from the 
menu and program execution begins with input of the system 


parameters. Input of the required matrices follows next. At 
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this point the user is prompted as to whether or not to have 
the modal distribution and gain matrices determined. Modal 
matrices can be used to determine the controlability and the 
observability/reconstructability of the system. Modal 


matrices are further explained in Section 7 to follow. 
b. Output Data Files 


The open-loop transfer function option outputs 
the system parameters, eigenvalues, transfer function gain, 
and order of the numerator to the OPTGRAPH file; OPTGROL 
data. This file may now be used in the OPTGRAPH FORTRAN 
program to obtain Pole-Zero and Root-Locus plots. Figures 
5.3 and 5.4 show OPTGRAPH's output for the open-loop 


transfer function. 
с. Solution Algorithms 


The eigenvalue analysis for the open-loop 
transfer function option is carried out using the same 
subroutines mentioned іп the discussion of the eigensystem 
analysis option. All eigen values are checked for system 
stability after being computed. If an eigenvalue falls 
outside the unit circle in the z plane, the system is 
unstable. In this case a statement of the system's insta- 
bility is sent to the screen and the transfer function 
analysis is bypassed. 

Transfer function analysis is carried out by 
OPTSYSX's subroutine TF. Option control flags ITFX and ITF2 
indicate to TF the type of transfer function and the anal- 
ysis desired. Subroutines POLES, ZEROS and RESID are called 
from TF to perform the analysis indicated by the flags. 


d. Example Problem Run 


The Kwakernaak and Sivan example problem is used 


to demonstrate the open-loop transfer function option. 
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TRANFUNC Example Run #2, Plot #1. 
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Figure 5.4 TRANFUNC Example Run #2, Plot #2. 
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Appendix F is а recording of a terminal session  demon- 
strating this option with the example problem data. Figures 
5.3 and 5.4 are the pole-zero and root-locus plots from 
OPTGRAPH for the example problem data. 


3. Noise Transfer Function Analysis 
a. Program Operation 


Figure 5.5 shows the Noise Transfer Function 
option menu. It includes the noise transfer function equa- 
tion and analysis options available. "Notes" are present in 
Figure 5.5 to remind the user of program restrictions and 
requirements germane to noise transfer function analysis. 

Program operation proceeds as in the open-loop 
case. Open-loop eigenvalues are calculated and the modal 
matrices are computed if desired. The closed-loop regulator 
dynamics matrix, (A-BF), is formed, and eigensystem analysis 
for the closed-loop follows. The closed-loop right eigen- 
vector matrix, (M) and closed-loop optimal regulator left 
eigenvector matrix, (M)”1, are computed. All eigenvalues are 
compared to the unit circle for system stability. 
Instability statements are issued if required and the 
transfer function analysis bypassed. Transfer function 
analysis is as discussed in the open-loop transfer function 


discussions. 
b. Output Data Files 


The noise transfer function option outputs the 
system parameters, eigenvalues, transfer function gain, and 
order of the numerator to the OPTGRAPH file; OPTGRNO data. 
This file may now be used in the OPTGRAPH FORTRAN program to 
obtain Pole-Zero and Root-Locus plots for the noise transfer 


function. 
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NOISE TRANSFER FUNCTION OPTIONS: 
Z/U = {H}*{ z{I}-{A}+{(B}*{F} }INV*{GAMD} 
OPTION 1 -- POLES, RESIDUES, AND ZEROS COMPUTED 
OPTION 2 -- ONLY POLES AND ZEROS COMPUTED 


1 
2 

OPTION 3 -- ONLY POLES AND RESIDUES COMPUTED 
4 


OPTION -- ЕХІТ -- RETURN ТО TRANFUNC OPTIONS 
SELECT OPTION 1l, 2, 3, OR 4 (SETS ITF2 = 1,2,3) 
NOTE: NOISE TF FUNCTION THRU CLOSED LOOP SYSTEM; 
NOTE: ODES AND ZEROS MUST BE SELECTED IF YOU 


WI TO USE OPTGRAPH. 
NOTE: REQUIRES НСО, NO>O, NG>O 





Figure 5.5 Noise Transfer Function Options. 


c. Solution Algorithms 


Analysis of the open-loop and closed-loop eigen- 
Systems analyzed in this option is done using the subrou- 
tines described in the open-loop transfer function 
discussion. Transfer function analysis is again carried out 
using the subroutine TF. This time the flags ITFX and ITF3 
are used to indicate the noise transfer function and desired 
analysis. Subroutines POLES,  ZEROS and RESID are called in 


TF to perform the desired analysis. 
d. Example Problem Run 


The Kwakernaak and Sivan example problem is used 
to demonstrate the noise transfer function option. Appendix 
G is а recording of а terminal session demonstrating this 


option with the example problem data. 
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4. Compensator Transfer Function Analysis 


a. Program Operation 


Figure 5.6 shows the Compensator Transfer 
Function option menu. It includes the compensator transfer 
function equation and analysis options available. "Notes" 


are present in Figure 5.6 to remind the user of program 
restrictions and requirements germane to compensator 
transfer function analysis. 

Program operation proceeds ina manner similar 
to the noise transfer function analysis. The open loop 
eigenvalues аге calculated and the modal matrices аге 
computed if desired. Then the closed-loop filter dynamics 


matrix, {A+KH} is formed.  Eigensystem analysis of the esti- 


mator follows. The closed-loop right eigenvector matrix, 
(M), the measurement  eigenvector, {H(BAR)*M}, and the 
closed-loop optimal filter left eigenvector matrix, {MY}, 


are computed. All eigenvalues are compared to the unit 
circle for system stability. Instability statements аге 
issued if required and the transfer function analysis 
bypassed. Transfer function analysis is carried out and the 


results printed as described in the the open-loop discus- 


sion. 
b. Output Data Files 
The compensator transfer function option outputs 
the system parameters, eigenvalues, transfer function gain, 


and order of the numerator to the OPTGRAPH file; OPTGRCM 
data. This file may now be used in the OPTGRAPH FORTRAN 
program to obtain Pole-Zero and Root-Locus plots for’ the 


compensator transfer function. 
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COMPENSATOR TRANSFER FUNCTION OPTIONS: 
ee ОО АСЕЕВ СЕО СНМ 
ОРТІОМ 1 -- POLES, RESIDUES, AND ZEROS COMPUTED 
OPTION 2 -- ONLY POLES AND ZEROS COMPUTED 


OPTION 3 -- ONLY POLES AND RESIDUES COMPUTED 
OPTION 4 -- EXIT -- RETURN TO TRANFUNC OPTIONS 
ОЕЕО РИО ТЕТО 3, OR 4 (SETSSITF3 = 1,2,3) 


NOTE: COMPENSATOR TRANSFER FUNCTION ABOVE IS FROM 
MEASUREMENT TO INPUT. 

NOTE: A COMPENSATOR TRANSFER FUNCTION MAY BE COMPUTED 
ONLY IF BOTH THE REGULATOR (F) AND FILTER (K) 
GAINS ARE AVAILABLE. THEY MAY BE CALCULATED IN 
THE "OPTREG" AND "KBFIL" PROGRAMS. 

NOTE: POLES AND pe. MUST BE SELECTED IF YOU WISH TO 
USE OPTGRAPH 

NOTE: REQUIRES NC»0, NO>0 


Figure 5.6 Compensator Transfer Function Options. 


c. Solution Algorithms 


Analysis of the open-loop and estimator eigen- 
systems in this option is done using the subroutines 
described previously discussed. Transfer function analysis 
is via the subroutine TF. This time the flags ITFX and ITF4 
are used to indicate the compensator transfer function and 
desired analysis. Subroutines POLES, ZEROS and RESID are 


called in TF to perform the desired analysis. 
d. Example Problem Run 


The Kwakernaak and Sivan example problem is used 
to demonstrate the compensator transfer function option. 
Appendix H is a recording of a terminal session  demon- 


strating this option with the example problem data. 
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5. Тһе MARKOV Parameter 


THIS OPTION DETERMINES THE CRITERIA FOR DECIDING WHEN 
A MARKOV PARAMETER IS ZERO. THE MARKOV PARAMETER 
INDICATES THE ORDER OF THE NUMERATOR POLYNOMIAL OF 
EACH TRANSFER FUNCTION. 


ALL "N" ZEROS OF THIS POLYNOMIAL ARE PRINTED OUT AND 
THIS TEST Rd .HOW MANY EXTRA ROOTS EXIST AT Z 0. 
LESS THAN 10.0*5*(-IE) IS CONSIDERED ZERO. 


THE DEFAULT VALUE OF THIS PARAMETER (IE) IS 6. 
IN OTHER WORDS, IE = 1.0E-6. 


IF YOU DESIRE A DIFFERENT MARKOV CRITERIA, 
ENTER THE INTEGER VALUE. 


IF YOU DESIRE THE DEFAULT VALUE, TYPE "6" 





Figure 5.7 MARKOV Parameter Option. 


6. НЕГР 


Figure 5.8 shows the HELP feature available іп 
TRANFUNC. The menu gives the user the option of obtaining a 
hardcopy of the information selected for viewing. 
Information available in HELP has been included in the 
terminal session for the eigensystem analysis. See Appendix 
Es 


7. Modal Matrices 


In all of the transfer function analysis options of 
TRANFUNC the user may select to have the modal matrices 


computed. Modal matrices can be used in determining whether 


a system is controllable and observable/reconstructable. I£ 
the modal control distribution matrix,  (TI*B) has no zero 
elements, the system is completely controllable. ТЕ the 
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"HELP" EIDEM OVERVIEW OF THE PROGRAM OPTIONS 
AVAILABLE IN  "TRANFUNC AND SOME OF THE EQUATIONS 
USED IN THE CALCULATIONS 


l -- TERMINAL VIEWING OF SELECTED INFORMATION 
- TERMINAL VIEWING PLUS HARDCOPY TO YOUR "A" DISK 


-- EXIT: RETURN TO THE TRANFUNC OPTIONS 
NOTE: OPTION "2" WILL PLACE THE FILE; "TRANFUNC HELP" 
ON YOUR А DISK CONTAINING THE ш» 
INFORMATION THAT YOU HAVE VIE 


ЕЕС OPTION: 1l, 2, OR 3 





Figure 5.8 TRANFUNC HELP Control Menu. 


modal measurement scaling matrix, (H(BAR)*T) has no zero 
elements, the system is completely observable/ 
reconstructable. With these matrices computed the user can 


tell at a glance whether the system is controllable апа 


observable with respect to any state/input. 
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VI. SPECIAL USE CAPABILITIES OF ORACLS 

The ORACLSX and TRANFUNC programs were designed with 
modular programing in mind. They have a main program which 
calls the selected analysis subroutine after system parame- 
ters and input data have been entered. This type of setup 
allows an interested user the freedom to develop an analysis 
subroutine Suited to a particular problem. This new subrou- 
tine could then be added to the main program as an addi- 
tional option. System parameters and data  input/output 
handling would still be handled by the main program. 
ORACLS's [Ref. 1] support and analysis subroutines could be 
used to develop this subroutine. This would minimize devel- 
opment time and any interface problems. Additional subrou- 
tines could be added as needed once interface problems were 
worked out. 

This capability is being used by а student working on a 
Masters Thesis in the area of self healing control systems 
at the Naval Postgraduate School, Monterey, California. It 
required a program which would input a large group of 
continuous matrices, discretize them, and then combine them 
into matrices appropriate for analysis in  ORACLS. The 
program was easily put together using the matrix handling 
subroutines of ORACLSX and ORACLS. Flexibility such as this 
gives ORACLS unlimited potential in the design and analysis 


of discrete control systems. 
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VII. CONCLUSIONS 


The advent of the digital computer and the micropro- 
cessor has moved the study and design of digital control 
systems to the forefront of control engineering. Control 
engineers and students are continually in need of computer 
programs to handle the analysis of discrete systems. The 
ORACLSX and TRANFUNC programs developed in this thesis give 
just such a capability to the CONTROLS analysis package. 
These programs when used in conjunction with the OPTSYSX 
program [Ref. 3], provide the user with the ability not only 
to handle analysis of continuous and discrete systems, but 
to be able to move from one to the other. Results from 
ORACLSX and TRANFUNC are compatible with the  OPTPLOT 
[Ref. 5] and OPTGRAPH [Ref. 4] programs. This allows a 
graphical displaying of the program results, a feature 
highly desirable in the study and design of controls 
systems. 


A few example problems were used to test program opera- 


tion and the accuracy of generated results. All tests 
produced favorable and accurate results. The user is 
cautioned, however that the programs developed in this 


thesis may not have been exercised for all cases of 


interest. While every effort has been made, within the time 
available, to ensure that the programs are free of computa- 
tional and logic errors, they cannot be considered vali- 


dated. Any application of these programs without additional 


verification is at the risk of the user. 
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APPENDIX А 
ORCONV EXAMPLE RUN 


BEGIN RECORDING OF TERMINAL SESSION 
рев 01/70.02 20:37: ]5 

controls 

DASD 200 LINKED R/O; R/W BY 0637P 

C (200) R/0 


PLOTTING IS DONE THRU DISSPLA. UP TO 20% OF YOUR DISK 
SPACE MAY BE NEEDED OTHERWISE ERROR MESSAGES 


wate als ate ate ale ate atc ats ats ale alc alc als atc atc alc ats ate at. Jes tes la wta ata nta nta nta nta nta ata a te ste ste’ ама ыб ы, > Ses Yo ate ste ste 2-2. ole als ate als ale ate ate ata m als ate ats ale 
0% 60% 6% 6% 9% 9% 69% 6% 06% 60% 6% 6% 6% 6% 6% 095 6% 2“ 0% 0% 6% 6% 6% 7 6% 4S 0S 4% 0% 4% оъ GX 9* 99 90V 4 2 0% 60% 6% 9% 9% $^ 


CONTROLS CONSISTS OF THREE INTERACTIVE PROGRAMS 


1 DACSAP SISO 
OPTSYS MIMO 
ORACLS MIMO 


HELP - PROGRAM DESCRIPTIONS 
SAMPLE PROBLEMS WITH PROGRAM DIRECTIONS 
E PM DOCUMENTATION SOURCES 


JONC UN 


ENTER |І, Ж». ма: ми: 4 5m E 


2% 2% 9% 9% 2% 9% 9% 9% 9% 9% 9% 9% 9% 2% 2% 9% 9% оъ 0S 0S ФО 2% 9% 0% оъ 2% 9% 9% 9% 9% 9% 2% 9% 9% 2% 9% 9% 9% 2% 9% 2% 9% 2% 9% оъ 9% 9% 9% 9% 9% 9% 9% оъ 9% 9% 2% 9% 9% 95 


NORMALLY IN ANY FORTRAN PROGRAM TWO NULL ENTREES TERMINATE 
THE PROGRAM SENDING YOU BACK TO THE CONTROLLING EXEC 


DESI pEPLACE R F (OS) NE 
(291) В 
"2517 REPLACES ' F (251) ' 
E251 ier 


-— = . е е е а «к, «аш, «жш, «ж, «к, еш, «нш; AO «шу «ы; «ен; «шу «ны; «шн; «ы; «к» «кы; «кыр GR «ны; GR «ы; тыз «с, тыз ет чыз “р т = а а е со =} =} «ы; «ы; «ы; «жш; «ж; «аш; «нш, «шн; «шн; «ж, е а а ао ено m» 


Jue THE ORACLS EXEC ALLOWS THE FOLLOWING OPTIONS: "eee 
1 ORACLSX FORTRAN: -------- DISCRETE SYSTEM ANALYSIS------- 
(A-D CONVERSION, OPT. REG., KALMAN-BUCY 
FILTER AND DISCRETE TRANSIENT ANALYSIS) 
2 TRANFUNC FORTRAN: -DISCRETE TRANSFER FUNCTION ANALYSIS- 
(SYSTEM EIGENVALUES, OPEN LOOP, NOISE 
AND COMPENSATOR T.F.s, MODAL MATRICES } 


3 ОРТР,ОТ FORTRAN: (PLOTTING OF TRANSIENT RESPONSES FOR 
STATES AND CONTROLS ) 


OPTGRAPH FORTRAN: ( POLE-ZERO, ROOT-LOCUS ) 

OUTPUT SELECTION -- TERMINAL OR "A" DISK FILE 

HELP -- PROGRAM AND DATA FILE RELATIONSHIP EXPLANATIONS 
EXIT --" KETURN ТО СМ5 


NSN сп л Lf 
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GENERAL ORACLS OPTIONS: 


1 -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
fee ORDSAL: DISCRETE SYSTEM ANALYSIS 
Bee SORTRAL: TRANSIENT ANALYSES OF A DISCRETE SYSTEM 
4 -- HELP: PROGRAM & EQUATION DESCRIPTIONS 
ee eee) RETURN TO ORACLS EXEC 
NOME: SELECT OPTION - 5; IF YOU ARE USING THE OUTPUT TO 


DISK OPTION AND би ULD LIKE TO CHANGE FILENAMES 


NOTE: NORMALLY ANY TIME DURING THE PROGRAM TWO SUCCESSIVE 


NULL "ENTERS" TERMINATES THE PROGRAM, SENDING YOU 
BACK TO THE ORACLS EXEC. 


SELECT AN OPTION: 1,2, 3, 4, ОВ 5. 


LIAE оса е е е е е е е е е е е е оо е се sw M M со е ом е е е ес он сом мо оо е со е е ee 


"HELP" PROVIDES, AN OVERVIEW OF THE PROGRAM OPTIONS 
AVAILABLE ІМ  "ORACLS" AND SOME OF THE EQUATIONS 
USED IN THE CALCULATIONS. 


1 -- TERMINAL VIEWING OF SELECTED INFORMATION 
2 -- TERMINAL VIEWING PLUS HARDCOPY TO YOUR "A" DISK 
3 -- EXIT: RETURN TO THE ORACLS OPTIONS 


NOTE: OPTION "2"eWELL PBACE THE FILE; "*ORACLS HELP" ON 
YOUR "A" DISK CONTAINING THE SELECTED INFORMATION 
THAT YOU HAVE VIEWED AT THE TERMINAL. 


БЕБЕСЕШӘРЛЕОЫЗ- 1, 2, ОК 3 


=» мис шо ао а ас а а еш» е е мо ее е е е еее е е е е е е е е е е е е е е TF Fe TCT е е е е ее е е ое е че ме ое е е е мо мо мм 


"HELP" INFORMATION OPTIONS: 

- OPTSYSX-ORACLS; EQUATION & MATRICES NOTATION COMPARED 
- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
- KBFIL: DISCRETE KALMAN-BUCY FILTER 
- OPTREG: DISCRETE OPTIMAL REGULATOR 
- ORTRAL: TRANSIENT ANALYSIS ОҒ А DISCRETE SYSTEM 
A EXIT: RETURN TO ORACLS OPTIONS 
- BRIT: RETURN TO ORACLS EXEC 

ШЕШЕ AN OPTION: 1,2; 3, 4, 5, 6, OR 7. 
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ORACLS IS A COMPLETELY INTERACTIVE DISCRETE SYSTEMS 
CONTROL PROGRAM. IT WILL SOLVE NUMEROUS CONTROL 
PROBLEMS OF THE FOLLOWING TYPES OF CONTROL EQUATIONS: 


CONTINUOUS SYSTEM:  (OPTSYSX NOTATION) 


SYSTEM EQUATION: XDOT = {F}*X +{G}*U +{GAM}*(W+WO) 

OUTPUT EQUATION: Y = (H)*X * (D)*U 

MEASUREMENT EQUATION: Z = {H}*X + (D}*U + V ; V:N(0,R) 

CONTROL LAW: U = -(C)%X 

COST EQUATION: J = 1/2*(INTGRL: {Y*{A}*Y + U*{B}*U} DT) 
DISCRETE SYSTEM: (ORACLS NOTATION) 


SYSTEM EQUATION: 
X((I*1)*T) = {A}*X(I*T) + (B}*U(I*T)+ (GAMD)*WD(I*T) 


OUTPUT EQUATION: Ү(І) = HFX) * (G)*U(I) 
MEASUREMENT EQUATION: Z(I) s (H)*X(I) * (G)*U(I) *V(I) 
CONTROL LAW: U = -{F}*X 
COST EQUATION: : X- х 

9 m Ју СЕН 1:175 SU s Ay 020 E 11; жыр CIS “Т) 


WHERE T IS THE SAMPLE TIME INTERVAL. 
"~a" = TRANSPOSE. 
WHERE "gen IS LIMIT AS "N" APPROACHES INFINITY. 


DO YOU WISH A LIST OF MATRICES USED IN OPTSYSX AND ORACLS? 
ANSWER (Y)ES OR (N) 


=» ае е е ес ме е е е е е е е е е е е е е е M M е е е е е мо е со е е соо ено е со е ено ме ом е се он «ы. «ы «ы» «ы «ж» «ы. «ы «ы «ж «шә «еш» «ы» «шы» 


MATRIX DESCRIPTION OPTSYSX NAME ORACLS NAME 


LIBE NE M M M е мо е е см е е е е е се е е е е е е е е е е е е -- се е е е е е е е е MU M чен алә ав ме мо мо мо ме мм ме 


SYSTEM MATRIX: dun 
CONTROL MATRIX: С 
NOISE MATRIX: ( GAM D) 
OBSERVABLES MATRIX: H 
MEASUREMENT-FEED FWD. DIST. MATRIX: D 

CONTROL GAIN MATRIX: C 

FILTER GAIN MATRIX: K 

OUTPUT COST MATRIX: А 

CONTROL COST MATRIX: B 

PROCESS NOISE MATRIX: 9 

MEASUREMENT NOISE MATRIX: 

CROSS PRODUCT WEIGHTING MATRIX: 
RICCATI EQUATION SOLUTION MATRIX: 
NOISE VECTOR: (И 
CONSTANT NOISE VECTOR: 


Ед СОК t> 
wy 


OO МАА 


"HELP" INFORMATION OPTIONS: 


l -- OPTSYSX-ORACLS; EQUATION & MATRICES NOTATION COMPARED 
2 -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
3 -- KBFIL: DISCRETE KALMAN-BUCY FILTER 
& -- OPTREG: DISCRETE OPTIMAL REGULATOR 
5 -- ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
ee EXIT: RETURN TO ORACLS OPTIONS 
EE SEXIT: RETURN TO ORACLS EXEC 
SELECT AN OPTION: 1,2, 3, 4, 5, 6; OR 7. 
2 


ORCONV IS AN ANALOG - DIGITAL MATRIX CONVERSION PROGRAM 
DANTEL ACCEPT INPUT FROM EITHER THE SCREEN OR FROM THE 
OPTMAT DATA" FILE FROM THE OPTSYSX PROGRAM. THE PROGRAM 
WILL CONVERT THE FOLLOWING CONTINUOUS SYSTEM MATRICES: 
ШІ 15), СВИ “tay, (Б), (Oy, (ТК) 
TO THE FOLLOWING DISCRETE SYSTEM MATRICES: 
CO БУУ СОЮ ЖО ТОКЕ Vl, (V2) 
AND ALSO COMPUTES THE DISCRETE CROSS PRODUCT WEIGHTING 
MATRIX; (S) 
WISH TO VIEW THE CONVERSION EQUATIONS USED IN ORCONV? 
ANSWER (Y)ES OR (N)O. 


THE FOLLOWING ARE THE CONVERSION EQUATIONS USED IN "ORCONV: 
ОА = xPCCF)*T) 


(B) = INTGRL:0-T ( (EXP({F}*TAU))*{G} ) DTAU 

{GAMD} = INTGRL:0-T ( EXP({F}*TAU)*{GAM} ) DTAU 

(Q) = INTGRL:0-T ( EXP({F}*TAU)*{A}*(EXP({F}-)*TAU) ) DTAU 
(В) » INTGRL:0O-T ( (B) * (H)-(TAU,O)*(A)*(H) (TAU,O) ) DTAU 


(Vl) » INTGRL:0-T ( EXP((F)-*TAU)*(Q)*(EXP((F))*TAU) ) DTAU 
{V2} = INTGRL:0-T ({R} + {H1}->(TAU,0)*{Q}*{H1}(TAU,0)) DTAU 


(S) s 2*INTGRL:0O-T ( (EXP((F)-*TAU))*(Q)*(H) (TAU,O) ) DTAU 
NOTE: + = | } TRANSPOSE 
Н}(ТТ, = INRGRL:0-TT ЕХР({Е ‘TAU G) ). DTAU 
ін Т ТТ 2) ZR INESRE 20 TT | КЕРЕГЕ ута 5 ) DTAU 
"TREIS TH SAMPLE TIME 


"HELP" INFORMATION OPTIONS: 

l -- OPTSYSX-ORACLS; EQUATION & MATRICES NOTATION COMPARED 
2 -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
3 -- KBFIL: DISCRETE KALMAN-BUCY FILTER 
4 -- OPTREG: DISCRETE OPTIMAL REGULATOR 
5 -- ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
6 -- EXIT: | RETURN TO ORACLS OPTIONS 
7 -- EXIT: RETURN TO ORACLS EXEC 

SELECT AN OPTION: 1,2, 3, 4, 5, 6, OR 7. 


6 
GENERAL ORACLS OPTIONS: 
1] -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
2 -- ORDSAL: DISCRETE SYSTEM ANALYSIS 
J -- ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
8 -- HELP: PROGRAM & EQUATION DESCRIPTIONS 
5 -- ЕЛІН RETURN TO ORACLS EXEC 
NOTE: SELECT ОРТТОМ - 5: ТЕ YOU ARE USING THE OUTPUT TO 


DISK OPTION AND WOULD LIKE TO CHANGE FILENAMES 

NOTE: NORMALLY ANY, TIME DURING THE PROGRAM TWO SUCCESSIVE 
NULL ENTERS  TERMINATES THE PROGRAM, SENDING YOU 
BACK TO THE ORACLS EXEC. 


SELECT AN OPTION: 1,2, 3, 4, OR 5. 


1 
ENTRR ТНЕ £ OF STATES (NS) OF THE SYSTEM MATRIX 
THE CONTINUOUS SYSTEM EQUATION: 
(М5) 
XDOT = (F)*X + (G)*U + (GAM)#(W+WO) 
(NS*NS) (NS*NC) (NS*NG) 
NS =? 
2 


HES HO THE КОП, SES CONTROLS (NC) OF THE CONTROL SYSTEM MODEL 


THE CONTINUOUS SYSTEM EQUATION: 
E 
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XDOT - (F)*X + (G)*U + (GAM)*(W+WO) 
{NS*NS} (NS*NC) (NS*NG) 


THE CONTINUOUS SYSTEM EQUATION: 
(NG) 
XDOT = {F}*X + {G}*U * (GAM)*(W*WO) 
(NS*NS) (NS*NC) (NS*NG) 


ENTE CIMA THE rix ОЕ LE oum OR OBSERVATIONS (NO) OF THE 


THE OUTPUT c - THE MEASUREMENT EQUATION: 
(NO) (NO) 
E AX + Dt) БЕ ооу 
(NO*NS) (NO*NC) (NO*NS) (NO*NC) 


ILL A FEED-FORWARD ae ee MATRIX 
{ D" - MATRIX} BE INPUT 


Gat OUTPUT EQUATION: THE MEASUREMENT EQUATION: 
Y = {H}*X + {D}*U ООО (руу кту 
TYPE YESKE TOR "NOS 
NOTE: SETS THE IFDFW FLAG 


WHAT IS THE SAMPLE TIME INTERVAL ("DELT") FOR GENERATING 
THE DISCRETE MATRICES? 


DELT = ? ENTER A DECIMAL ANSWER. 
ШЕСІ OR S5 OR 05 ETC) 


FLAG/PARAMETER SETTINGS FOR THIS RUN ARE AS FOLLOWS: 


F G GAM À B 8 К -- (CONTINUOUS NAME ) 
І 1 0 0 О О О -—- (MATRICES NAME) Ер) 
А В САМР 0 R Ve Zee JO DISCRETE NAME 
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SAMPLE TIME = О и 
ORDER OF SYSTEM = 2 


NUMBER OF CONTROLS = 1 
NUMBER OF OBSERVATIONS = 0 
NUMBER OF PROCESS NOISE SOURCES = 0 


== = ee SS eee -а а а  - а нт «ж нт «ж |. — |. ы ан сы» «шқ 


ENTER THE SYSTEM MATRIX {"'F"-MATRIX} 
DIMENSION = # STATES {NS} X # STATES {NS} 
YOU MAY ENTER THE MATRIX OR HAVE THE PROGRAM GENERATE 
A DESIRED MATRIX BY SELECTING ONE OF THE FOLLOWING: 
1 - NULL MATRIX 


2 - IDENTITY MATRIX 
3 - IDENTITY MATRIX SCALED BY A DESIRED VALUE 
4 - ANY OTHER MATRIX 
SELECT ANY OPTION 
4 
y THE ELEMENT F( 1, 1)= 
THE ELEMENT F( l, 2)= 
o THE ELEMENT IU. je 
_„ g THE ELEMENT Е( 2, 2)= 
THE SYSTEM MATRIX {"F'-MATRIX}... 
0.00000 1.00000 
0:00000 -4.60000 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
А TYPE "YES" OR "NO". 


LIE ONE а «ж, «жа «жа, «ж, o е е е е е е е е е ес е а оо е е чы ==» з ч чыз == жы тз чы ч ч ч ч == се 7^ 171 cl M M M M M X M M ч ч оо оо е ть 


ENTER THE CONTROL DISTRIBUTION MATRIX {"G"-MATRIX}. 
DIMENSION = # STATES {NS} X # CONTROLS {NC} 
YOU MAY ENTER THE MATRIX OR HAVE THE PROGRAM GENERATE 
A NULL MATRIX BY SELECTING ONE OF THE FOLLOWING: 
1 - NULL MATRIX 
2 - ANY OTHER MATRIX 
SELECT AN OPTION 


2 


= æ мә ч. мь  - - м | жы а 
е аә сә еә ољ еә еә аә еә еә еә мо е е Әәә е “е. з з ч е «ез «е» ғ» ез әз с» ә.з s“ чә ы .ь -ь “< ев шә =з =з шэ шэ шэ = 
=» -- еә w м = 


THE ELEMENT G( 1, 1)= 
THE ELEMENT G( 2, 1)- 


787 
THE CONTROL DISTRIBUTION MATRIX ('"G'"-MATRIX) 
0.00000 
0.78700 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
ТҮРЕ YES LOR “МО! 


= --. =» =» мо ме е мо мо е ме - е е мо е мо мо мо не мо мо мо ме мо o o зо -. -. ы. ы» ы. но ә. -. мо -. ӘӘ Әһ Әһ ы» о -о ыз ма ы» ха -- - eel el ме 


рО YOU WISH TO GENERATE AN "OPTMAT' DATA FILE TO SAVE THE: 


ШІ (5Ү65ТЕМ), (G) (CONTROL), 
{H} (OBSERVABLES), | (D) (MEASUREMENT- FEED FWD. DIST.), 
(GAM) (NOISE), (A) (OUTPUT COST), 


(B) (CONTROL COST), (Q) (PROCESS NOISE INTENSITY (PSD)), 
(R) (MEASUREMENT NOISE INTENSITY (PSD)) 


MATRICES AS YOU HAVE JUST a FOR REENTRY TO THE 
OPTSYSX OR ORACLS PROGRAM 


HOTE PHIS WILE WRITE OVER THE OLD "OPTMAT' FILE IF YOU 
HAVE USED ONE FOR THIS RUN. 


(Y OR N) 


= æ ми ес ел е ел е е е е е се е е е е ғ. әль е е се оо е со се «ә» мо сы» ар а е фев: со сыз сыз сыз» е е сыз сыз мо ее е се е се се 1 сәз е ыз: OM M M 27 


ОРТМАТ DATA FILE INPUT TO ORCONV PROGRAM: ORACLES PROGRAMS 


F MATRIX 2 ROWS 2 COLUMNS 
0.0000000E+00 1.0000000E+00 
0.0000000E+00 -4.6000000Е%00 

G MATRIX 2 ROWS 1 COLUMNS 
0.0000000E+00 
7.8700000Е-01 


-- мо мо ме е ме ме мо мм е ме е е е е е е а е е ә» чи» е.» чи» сыз оо ә е с.» е сы» сә» ыз сыз сыз сыз сыз оо мн е е е сыз мо ыз ма ә е е е е -. е е ми 


ORCONV: ANALOG - DIGITAL CONVERSION...DISCRETE MATRICES... 


THE SAMPLE TIME INTERVAL FOR THIS RUN - 0.1000 
А МАТКІХ 2 ROWS 2 COLUMNS 
1. 0000000Е+00 015572 Ë 02 
0.0000000E+00 Е ЭЕ 
В МАТКІХ 2 ROWS l COLUMNS 
DE ES0U239E-03 
595097559 E-02 
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-. -.а ме ме е е е е е е е - мо е а мо  - 5 ар -Ó е е е е се е е се се е е е е се се ее се е е е ә. -. - 5 ә. ә мн е е мә. к. ә. әс ы. ж = 


DO YOU WISH TO GENERATE AN OPTMATD DATA FILE ТО 5АУЕ THE: 


(A) (SYSTEM), (B) (CONTROL), 
(H) (OBSERVABLES), (G) K(MEASUREMENT- FEED FWD. DIST.), 
(GAMD) (NOISE), (F) (CONTROL GAIN), 


(К) (FILTER CAIN), СО, О ОШЕН ШЫ, 

(R) (CONTROL COST), (V1) (PROCESS NOISE INTENSITY (PSD)), 
(V2) (MEASUREMENT NOISE INTENSITY (PSD)), 

(S) (DISCRETE CROSS PRODUCT WEIGHTING) 

MATRICES FOR REENTRY 19 ІНЕ 0РТ5Ү5Х OR ORACLS PROGRAMS? 


GENERAL ORACLS OPTIONS: 
1 -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
2 -- ORDSAL: DISCRETE SYSTEM ANALYSIS 
J -- ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
4 -- ПЕРЕН PROGRAM & EQUATION DESCRIPTIONS 
Э ти: RETURN TO ORACLS EXEC 
NOTE: SELECT OPTION - 5; IF YOU ARE USING THE OUTPUT TO 
DISK OPTION AND WOULD LIKE TO CHANGE FILENAMES 
NOTE: NORMALLY ANY TIME DURING THE PROGRAM TWO SUCCESSIVE 
NULL "ENTERS" TERMINATES THE PROGRAM, SENDING YOU 
BACK TO THE ORACLS EXEC. 


SELECT AN OPTION: 1,22 9. ЖҮРЕ. 


...-ORACLS IS NOW TERMINATED. ...RETURNING TO ORACLS EXEC... 


MESSAGE SUMMARY: MESSAGE NUMBER - COUNT 
187 22 
елу THE ORACLS EXEC ALLOWS THE FOLLOWING OPTIONS: ЖЛ 
1 ОРАС &Х КОКТЕАМ ш DISCRETE ЗУФШЕМ АМЛИМВНС5- oa 
(A-D CONVERSION . REG.. KALMAN-BUCY 
FILTER AND DISCRETE TRANSIENT ANALYSIS) 
2  TRANFUNC FORTRAN: -DISCRETE TRANSFER FUNCTION ANALYSIS- 
(SYSTEM EIGENVALUES, OPEN LOOP, NOISE 
AND COMPENSATOR T.F:s, MODAL MATRICES} 


3 OPTPLOT FORTRAN: (PLOTTING OF TRANSIENT RESPONSES FOR 
STATES AND CONTROLS ) 


4 OPTGRAPH FORTRAN: ( POLE-ZERO, ROOT-LOCUS ) 
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5 OUTPUT SELECTION -- TERMINAL OR "A" DISK FILE 

6 HELP -- PROGRAM AND DATA FILE RELATIONSHIP EXPLANATIONS 
КОШЕ et -- RETURN TO CMS 

ж SELECT AN OPTION, ANY OTHER INPUT RETURNS YOU TO MENU ** 


ж = — eee eee eee 


record off 
END RECORDING OF TERMINAL SESSION 
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APPENDIX B 
KBFIL EXAMPLE RUN 


==» еә — и. əм -. ә.  - ғ. -. -. - -ы, «ме ың -- ңын «м ә аә ә әәә ә ә Б Б ә — =— ӘБ «ж» Әә о ә ӘБ әәә ә «м. әже -ы, ән, ән, оғ, GO RD ғи ән ә il m 


"HELP" INFORMATION OPTIONS: 


L -- OPTSYSX-ORACLSX; EQUATION & MATRICES NOTATION COMPARED 
2 -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
3 -- KBFIL: DISCRETE KALMAN-BUCY FILTER 
4 -- OPTREG: DISCRETE OPTIMAL REGULATOR 
5 -- ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
6 -- EXIT: RETURN TO ORACLSXK OPTIONS 
7 -- EXIT: RETURN TO ORACLS EXEC 
SELECT AN OPTION: 1,2, 3, 4, 5, 6, OR 7. 
3 


KBFIL WILL SOLVE THE DISCRETE TIME-INVARIANT OPTIMAL 
KALMAN-BUCY FILTER PROBLEM. SOLUTION OF THIS PROBLEM 
PRODUCES THE OPTIMAL FILTER GAIN, "K" MATRIX WHICH MAY 
BE OUTPUT TO THE OPTMATD FILE FOR USE IN THE TRANFUNC 
PROGRAM THE OPTIMAL ESTIMATOR PROBLEM IS TO CONSTRUCT 

AN ESTIMATE, XHAT(I) OF X(I) FROM KNOWLEDGE OF THE 
OUTPUTS, Y(I) SUCH THAT THE COST CONSTRAINT IS MINIMIZED. 
IF THE SYSTEM IS OBSERVABLE/RECONSTRUCTIBLE AND 
CONTROLLABLE THE ASYMPTOTIC OPTIMAL OBSERVER PROBLEM 
EXISTS AND XHAT(I) IS GIVEN BY: 


XHAT (I*1)-(A) *XHAT(I)* (B)JU(I)* (K) * (Y(I)- (H) *XHAT(I)- (G)U(I)) 
WITH FILTER GAIN 

(K) = MAJAP] (H) (R) > ТЕРЕ ЕНЕ 
AND WITH (P) SATISFYING: 


(Ploss(M)*CELZ(M)- + (КҮЗЕТТЕ Re ry O (Con). 
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FOR 
СЕ ЕСЕ) 
THE MATRIX "P" REPRESENTS THE RECONSTRUCTION ERROR. 


WHERE: ( )-^ = MATRIX TRANSPOSE 
"HELP" INFORMATION OPTIONS: 
l -- OPTSYSX-ORACLSX; EQUATION & MATRICES NOTATION COMPARED 
2 -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
3 -- KBFIL: DISCRETE KALMAN-BUCY FILTER 
4 -- OPTREG: DISCRETE OPTIMAL REGULATOR 
25 -- ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
6 -- EXIT: | RETURN TO ORACLSX OPTIONS 
7 -- EXIT: RETURN TO ORACLS EXEC 
SELECT AN OPTION: 1,2, 3, 4, 5, 6, OR 7. 


6 
GENERAL ORACLSX OPTIONS: 
1 -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
Zea ORDSAL: DISCRETE SYSTEM ANALYSIS 
3 -- ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
4 -- HELP: PROGRAM & EQUATION DESCRIPTIONS 
эу „== ол RETURN TO ORACLS EXEC 
SELECT AN OPTTON: 1.2. 3, 4, ОК 5. 
NOTE: SELECT OPTION - 5; a YOU ARE USING THE OUTPUT TO 


DISK OPTION AND POUL LIKE TO CHANGE FILENAMES 

NOTE: NORMALLY ANY. ІГГІЛЕ DURING THE PROGRAM TWO SUCCESSIVE 
NULL "ENTERS" TERMINATES THE PROGRAM, SENDING YOU 
BACK TO THE ORACLS EXEC. 


ORDSAL OPTIONS: 
-- KBFIL : DISCRETE KALMAN-BUCY FILTER ANALYSIS 
-- OPTREG: DISCRETE OPTIMAL REGULATOR ANALYSIS 
EXIT: RETURN TO ORACLSX OPTIONS 
== (ЕМЕ RETURN TO ORACLS EXEC 
SELECT AN OPTIONER, 2, 3, ORA 


£ € hN n 
! 
| 
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NOTE: TO PERFORM COMPENSATOR TRANSFER FUNCTION ANALYSIS ІМ 

TRANFUNC REQUIRES BOTH THE {K} AND (F) MATRICES. DO 

OPTREG" FIRST AND THEN "КВЕТ AS THIS ENSURES THE 
CORRECT (H) MATRIX WILL BE AVAILABLE IN TRANFUNC. 


THE ҒА 1 "B" "САМР" "Hu "G "Ez K t " DE "үк м2" AND еп 


matrices FROM YOUR’ PREVIOUS ORACLS RUN WERE SAVED. 
THE FOLLOWING OPTIONS ARE AVAILABLE: 
1. USE ALL OF THE SAME MATRICES AGAIN. 
2. USE SELECTED MATRICES AGAIN. 
3. INPUT ALL NEW MATRICES. 

ENTER 1, 2, OR 3. 
NOTE: EACH SAVED MATRIX WILL BE REDISPLAYED AT 

THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 


NOTE: A CHANGE IN THE ORDER OF THE SYSTEM; -NS- 
REQUIRES SELECTION OF OPTION 5 


DO YOU WISH TO SAVE THE "A''-MATRIX FROM THE LAST 
RUN TO BE USED IN THE FOLLOWING RUN? 


NOTE: THE MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 

TYPE "YES" OR "NO". 


ENTER_THE f OF PROCESS NOISE SOURCES (NG) OF THE 
("GAMD'"-MATRIX). 


THE DISCRETE SYSTEM EQUATION: 
(NG) 
X(I*1) » (A)*X(I) * (B)*U(I) + (GAMD)** (WD* WO) 
(NS*NS) (NS*NC) (NS*NG) 


ENTER THE OF MEASUREMENTS OR OBSERVATIONS (NO) OF THE 
{( H -MATRIX 


THE OUTPUT EQUATION: THE MEASUREMENT EQUATION: 
(NO) (NO) 
YLI) = (H) X(T) (G) US zl) - (H)*X(I) * (G)*U(I) * V(I) 
(NO*NS) (NO*NC) (NO*NS) (NO*NC) (NO*1) 
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2 
THE SYSTEM MATRIX ("А"”-МАТЕШХ)... 
1.00000 0.08016 
0.00000 0.63128 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
ТҮРЕ YES OR "NOIR 
n 


тә ML ML. о» ао шы с» о: а: ее ов ә» ә» е» == с» «са о» о» он о: өң е» —————— 


ENTER THE PROCESS NOISE DISTRIBUTION MATRIX ("GAMD'-MATRIX). 
DIMENSION = # STATES {NS} X # PROCESS NOISE SOURCES {NG}. 


YOU MAY ENTER THE MATRIX OR HAVE THE PROGRAM GENERATE 
A DESIRED MATRIX BY SELECTING ONE OF THE FOLLOWING: 


1 - NULL MATRIX 
- IDENTITY MATRIX 
- IDENTITY MATRIX SCALED BY A DESIRED VALUE 
- ANY OTHER MATRIX 
SELECT ANY OPTION 


£ Cc N 


= = = но е е е е е е == е е е е е е ето ое е е е е е е е е е е е е ма ее чо е е е ае е е еше еш» е е е е е е е е ом е ое е ма ом 


THE PROCESS NOISE DISTRIBUTION MATRIX ("GAMD'-MATRIX)... 


1.00000 0.00000 
0.00000 1.00000 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
ТҮРЕ YES OR UNO”, 


n 
ENTER THE MEASUREMENT SCALING MATRIX ("H'-MATRIX). 
DIMENSION = # OBSERVATIONS {NO} X # STATES (NS). 
YOU MAY ENTER THE MATRIX OR HAVE THE PROGRAM GENERATE 
A DESIRED MATRIX BY SELECTING ONE OF THE FOLLOWING: 
1 - NULL MATRIX 
2 - IDENTITY MATRIX 
3 - IDENTITY MATRIX SCALED BY А DESIRED VALUE 
& - ANY OTHER MATRIX 
SELECT ANY OPTION 
2 


— — < == = ао мо а а ма ма ме ма шема еее е е мо е е е е е есен а е ее е е е е е ее е е е се е е е соо мес е ос е е е ен е р 


THE MEASUREMENT SCALING MATRIX ("H'-MATRIX). 


1.00000 0.00000 
0.00000 1.00000 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" ORP NOD P 


Со е мо ес е еее е е ме е е е е ом а е е е е е е е е е е е е е ее е се е е е е е е е е е е е е е ее е е р әр әәә ел» еше өше нә 


ENTER THE PROCESS NOISE INTENSITY (PSD) MATRIX 
ССМ ИЕ а 


DIMENSION = # PROCESS NOISE SOURCES (NG) X 
# PROCESS NOISE SOURCES (NG). 


YOU MAY ENTER THE MATRIX OR HAVE THE PROGRAM GENERATE 
A DESIRED MATRIX BY SELECTING ONE OF THE FOLLOWING: 


1 - NULL MATRIX 
2 - IDENTITY MATRIX 
3 - IDENTITY MATRIX SCALED BY A DESIRED VALUE 
4 - ANY OTHER MATRIX 
SELECT ANY OPTION 
2 
THE PROCESS NOISE INTENSITY (BSD) MATRIX.. m.. Vl 
80898 08088 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE YES OR КО 


атр мо а ес е е е е е е е е е е е е се е со «к» «к» «к» е е е е е е е е е е е е е е е е е е е е ое сее е е се е е е е со е е еке еше өш» 


ЕМТЕК ТНЕ MEASURED S Y НІ ЕТ раа (PSD) MATRIX 


DIMENSION = # OBSERVATIONS (NO) X # OBSERVATIONS (NO) 
Á DESIRED MATRIX BY SELECTING ONE OF THE FOLLOWING? ы 
1 - NULL MATRIX 
2 - IDENTITY MATRIX 
3 - IDENTITY MATRIX SCALED BY A DESIRED VALUE 
4 - ANY OTHER MATRIX 
SELECT ANY OPTION 


THE MEASUREMENT NOISE INTENSITY (PSD) MATRIX....V2.... 
1.00000 0.00000 
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0.00000 1.00000 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
ТҮРЕ ЧАБЫ ШОР "О". 


LEE E E = ë = -. -- — = =» =» аә» == aoe æ =» = 
=» =» = m ccm еә ср e a алә e ë 10" LS Ne =» =» еә а е е е - ав ми 


FLAG/PARAMETER SETTINGS FOR THIS RUN ARE AS FOLLOWS: 


THE SAMPLE TIME INTERVAL FOR THE SYSTEM IS = 0.10000 
ORDER OF SYSTEM = 2 
NUMBER OF CONTROLS = 1 
NUMBER OF OBSERVATIONS = 2 
NUMBER OF PROCESS NOISE SOURCES = 2 
СОЈА ЕВОСҮУ КІБТЕК АМАПГ,Ү515................ ORACLES PROGRAMS 


S SAMPLE TIME INTERVAL (DELT) FOR THIS SYSTEM IS = 0.100 


ASYMFI: SUBROUTINE TO SOLVE THE DISCRETE INFINITE-DURATION 
OPTIMAL FILTER PROBLEM. 


A MATRIX 2 ROWS 2 COLUMNS 
1. 0000000Е+00 S085» 72 91 02 
0.0000000E-*00 2517930651 208 
САМО MATRIX 2 ROWS 2 COLUMNS 
1. 0000000Е+00 0000000Е+00 
0.0000000Е%00 1.0000000Е+00 
Н МАТКІХ 2 ROWS 2 COLUMNS 
1. 0000000Е+00 0.0000000Е%00 
О. 0000000Е+00 1. 0000000Е+00 
INTENSITY MATRIX FOR COVARIANCE OF MEASUREMENT NOISE 
V2 MATRIX 2 ROWS 2 COLUMNS 
1.0000000E-*00 0.0000000Е%00 
0.0000000Е»00 1.0000000Е%00 
Vl MATRIX 2 ROWS 2 COLUMNS 
. 0000000Е+00 0. 0000000Е+00 
0.0000000Е%00 1.0000000Е+00 
INTENSITY MATRIX FOR,COVARIANCE OF PROCESS NOISE; 
{GAMD}*{V1}*{GAMD} 
MATRIX 2 ROWS 2 COLUMNS 
1.0000000Е»00 0.0000000Е%00 
0. 0000000Е+00 1.0000000Е+00 
KALMAN-BUCY FILTER GAIN 
K MATRIX 2 ROWS 2 COLUMNS 
6 1 0913667E-01 t .9392360E 02 
5 1]7731EÉE-03 3.4673134Е-01 


STEADY-STATE VARIANCE MATRIX OF RECONSTRUCTION ERROR 
Р МАТКІХ 2 ROWS 2 COLUMNS 
1-762 30950 E: 00 3.1174278E-02 
3.1174278Е-02 1.2188858E+00 
EIGENVALUES OF P 
EVLP | MATRIX 2 ROWS 1 COLUMNS 
: 2164957E+00 
.6254851Е+00 
A-KH MATRIX 2 ROWS 2 COLUMNS 
3.8086333E-01 3.0775569EÉE 0° 
-3 590177 1E Oj 2.8&55231F 01 
EIGENVALUES OF A-KH MATRIX 
2.8564526Е-01 0.0000000Е%00 
7977038Е-01 0.0000000Е%00 


DO YOU 
{А} 
(H) 
{САМР} 
(K) 
СК) 
v2} 


WISH TO GENERATE AN "OPTMATD" DATA FILE TO SAVE THE: 
(SYSTEM), {B} (CONTROL), 

(OBSERVABLES), {С} (MEASUREMENT- FEED FWD. DIST.), 
(NOISE), (F) (CONTROL GAIN), 


(FILTER GAIN)  (Q) (OUTPUT COST), 
(CONTROL COST) (V1) (PROCESS NOISE INTENSITY (PSD)) 
(MEASUREMENT NOISE INTENSITY (PSD)) 


(S) (DISCRETE CROSS PRODUCT WEIGHTING) 
MATRICES FOR REENTRY TO THE OPTSYSX OR ORACLS PROGRAMS? 


Tol. 
1 - 
2 - 
3 - 
a 
КОТЕ: 
4 


(Y OR N) 


еа сна «на «жау «жа «жа «жа «ж «ла «жа «на «жа «ша «ша «ша «кш. «ша оа е е е мо ес е «ик: «жа «ша «жа «жа «ша «жа «жа «к. «жа е е е е е е е е е со 1 WW M 1 е е е «ша 


ORDSAL OPTIONS: 
- KBFIL : DISCRETE KALMAN-BUCY FILTER ANALYSIS 
- OPTREG: DISCRETE OPTIMAL REGULATOR ANALYSIS 
- EXIWE RETURN TO ORACLSX OPTIONS 
= JEXT E: RETURN TO ORACLS EXEC 
SELECT ANĦOPTION: 1, 2. 2. ORTI. 
RARFSRUNIGOYRENSUTOR TANGER, RUNCETON, AMMETSTS н 


"OPTREG" FIKST AND THEN "KBFIL AS THIS ENSURES THE 
CORRECT (H) MATRIX WILL BE P AVAILABLE IN TRANFUNC. 


- а ме ое е е е е е е е е е см е е со е е е е е е е е е е е ес мн ес е ес е = «жа «жа «ша «ша «ла «ла «ша «нша «ла оо се со е е е е еа 


..ORACLS IS NOW TERMINATED..... RETURNING TO ORACLS EXEC... 
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APPENDIX C 
OPTREG EXAMPLE RUN 


= o— жо. — — — — ғ. -а шю шю чє m эъ чш == шы е е е е см е е е е е е е е е е ее е «же «ж «ше шы == -. .-. тє шы он» = ==  — = © ыа шю шю эы = мм 


"HELP" INFORMATION OPTIONS: 


l -- OPTSYSX-ORACLS; EQUATION & MATRICES NOTATION COMPARED 
2 -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
3 -- KBFIL: DISCRETE KALMAN-BUCY FILTER 
4 -- OPTREG: DISCRETE OPTIMAL REGULATOR 
5 -- ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
6 -- EXIT: | RETURN TO ORACLS OPTIONS 
7 -- EXIT: RETURN TO ORACLS EXEC 
SELECT AN OPTION: 1,2, 3, 4 Ps pr 
4 


OPTREG IS A DISCRETE SYSTEM ANALYSIS PROGRAM THAT WILL 
ACCEPT INPUT FROM EITHER THE SCREEN OR FROM AN "OPTMATD" 
DATA FILE. IT USES THE DISCRETE MATRICES: (A), (B), (H), 
(Q), AND (R) TO COMPUTE THE OPTIMAL FILTER GAINS; (FHAT) 
AND (F). IF THE SYSTEM HAS A CROSS PRODUCT WEIGHTING MATRIX 
(S), THE PROGRAM WILL CALCULATE A NEW GAIN MATRIX; (FE) 
WHICH WILL ELIMINATE THE CROSS PRODUCT TERM; (S) IN THE 
QUADRATIC SCALAR FUNCTION. NEW MATRICES; (AHAT), (BHAT), 
(QHAT), AND (RHAT) ARE THEN CALCULATED AND USED AS INPUTS 
TO THE PROGRAM. 

DO YOU HAVE A CROSS PRODUCT WEIGHTING MATRIX ? 

ANSWER (Y)ES OR (N)O. 
"HELP" INFORMATION OPTIONS: 

l -- OPTSYSX-ORACLS; EQUATION & MATRICES NOTATION COMPARED 
2 -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
3 -- KBFIL: DISCRETE KALMAN-BUCY FILTER 
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-- OPTREG: DISCRETE OPTIMAL REGULATOR 
ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
-- EXIT: | RETURN TO ORACLS OPTIONS 
-- EXIT: | RETURN TO ORACLS EXEC 
SELECT AN OPTION: 1,2, 3, 4, 5, 6, OR 7. 


ч боол Es 
! 
! 


6 
GENERAL ORACLS OPTIONS: 
1 -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
2 -- ORDSAL: DISCRETE SYSTEM ANALYSIS 
3 -- ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
& -- НЕГР: PROGRAM & EQUATION DESCRIPTIONS 
EN ЕЛШІ RETURN TO ORACLS EXEC 
NOTE: SELECT OPTION - 5: IF YOU ARE USING THE OUTPUT TO 


DISK OPTION AND WOULD LIKE TO CHANGE FILENAMES 

NOTE: NORMALLY ANY, TIME DURING THE PROGRAM TWO SUCCESSIVE 
NULL "ENTERS" TERMINATES THE PROGRAM, SENDING YOU 
BACK TO THE ORACLS EXEC. 


SELECT AN OPTION: 1,2, 3, 4, ОВ 5. 


-. ә -ә -ә ме оа е е е оо е е е зы» зы» әр «к. зы» ек» еш» ек» ек» «кш» «к» «ше» е со мо мо ән ә ән е е е е а е е е е е е е со ое е е оо е е е е е «шә шо оо е 


ORDSAL OPTIONS: 

-- KBFIL : DISCRETE KALMAN-BUCY FILTER ANALYSIS 
-- OPTREG: DISCRETE OPTIMAL REGULATOR ANALYSIS 
EXIT: RETURN TO ORACLS OPTIONS 

== ҮҮТ: RETURN TO ORACLS EXEC 

Seber АН ОРТТОК“ 1. 2. 3. OR 4- 

NOTE: TO PERFORM COMPENSATOR TRANSFER НЕ ӨН ANALYSIS ІМ 

TRANFUNC REQUIRES BOTH THE {K}, AND {F} MATRICES. DO 


"ОРТКЕС' FIRST AND THEN "KBFIL" AS THIS ENSURES THE 
CORRECT (H) MATRIX WILL BE AVAILABLE IN TRANFUNC. 


T Ç кн 
! 
! 


аә — = = єз ©з сш тыз мл е е е ее е е е е е е е е ме се е е е жыз ты. == == жы} е е е ос е чы ты» чыз чыз «кә чыъ =з жыз чы» =з =з ©з ч со е == е ое е оз 


ТН EA TAN ГСАМр” "н", O tt Е" "к " " р r VE; Ve AND 
n MATRICES FROM YOUR’ PREVIOUS BRACES RU RUN WERE SAVED. 


THE FOLLOWING OPTIONS ARE AVAILABLE: 
1. USE ALL OF THE SAME MATRICES AGAIN. 
2. USE SELECTED MATRICES AGAIN. 
3. INPUT ALL NEW MATRICES. 


ENFER 1, 2; ОК 3. 
NOTE: EACH SAVED MATRIX WILL BE REDISPLAYED AT 


THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
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INDIVIDUAL MATRIX ELEMENTS. 


NOTE: A CHANGE IN THE ORDER OF E SYSTEME -NSS 
REQUIRES OPTION 3 SELECTION 


DO YOU WISH TO SAVE THE "A"-MATRIX FROM THE LAST 
RUN TO BE USED IN THE FOLLOWING RUN? 


NOTE: THE MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SE UENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS ` 

TYPE "YES" OR "NO". 


DO YOU WISH TO SAVE THE "B'-MATRIX FROM THE LAST 
RUN TO BE USED IN THE FOLLOWING RUN? 


NOTE: THE MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 

TYPE "YES" OR "NO". 


NIER THE if OF PROCESS NOISE SOURCES {NG} OF THE 
GAMD -MATRIX}. 


THE DISCRETE SYSTEM EQUATION: 
(NG) 
X(I+1) = {A}*X(I) + (B}*U(I) + {(GAMD}*(WD+WO) 
{NS*NS) {NS*NC} (NS*NG) 


ep] 


СМЕНЕТИ THE EOM OF MEASUREMENTS OR OBSERVATIONS (NO) OF THE 


THE OUTPUT EQUATION: THE MEASUREMENT EQUATION: 

(NO) (NO) 

ҮСЕ СОЮ Z[1)=(H)#K(I) СЕЕ ОӨЖ (Тој 
(NO*NS) (NO*NC) (NO*NS) (NO*NC) (NO*1) 


-— сə» чо ее а е а е е е е е е е е е е е е е е е оо е еа е мо мо со со TH ww ws е е мо мо е е е се е е а е BW а е ме оо со со се се е ее сее се 


IS THE SYSTEM YOU WISH TO EVALUATE A STABLE SYSTEM ? 
(Y)ES OR (N)O ? 
NOTE: ANSWER "NO" IF IT IS UNSTABLE OR YOU ARE NOT SURE. 
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"OPTREG" WILL TEST THE MATRIX (A-BF) FOR STABILITY 
RELATIVE TO .9. IF A STABILIZING GAIN IS a E 
IT WILL BE COMPUTED AND NOTED IN THE RESULTS. 


`... sss 
THE SYSTEM MATRIX ( А"-МАТЕТХ)... 
E0980 — 20904 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
E ТУРЕ YES OR "МО". 
THE CONTROL DISTRIBUTION MATRIX ("B'-MATRIX)... 
89349 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
- INDE C YES" OR "NO". 
ENTER THE MEASUREMENT SCALING MATRIX ("H'-MATRIX). 
DIMENSION = # OBSERVATIONS (NO) X # STATES (NS). 
A R O 
1 - NULL MATRIX 
2 - IDENTITY MATRIX 
3 - IDENTITY MATRIX SCALED BY A DESIRED VALUE 
4 - ANY OTHER MATRIX 
SELECT ANY OPTION 
2 
THE MEASUREMENT SCALING MATRIX ("H'"-MATRIX)... 
09098 — 008099 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
Ë PFPE YES" OR “NO”. 


ENTER THE OUTPUT MEASUREMENT COST MATRIX ("Q'-MATRIX). 
DIMENSION = # OBSERVATIONS {NO} X # OBSERVATIONS {NO}, 


YOU MAY ENTER THE MATRIX OR HAVE THE PROGRAM GENERATE 
A DESIRED MATRIX BY SELECTING ONE OF THE FOLLOWING: 


1 - NULL MATRIX 
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2 - IDENTITY MATRIX 
3 - IDENTITY MATRIX SCALED BY A DESIRED VALUE 
4 - ANY OTHER MATRIX 

SELECT ANY OPTION 


2 
THE OUTPUT MEASUREMENT COST MATRIX {"Q"-MATRIX}... 
1.00000 0.00000 
0.00000 1.00000 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
ТҮРЕ YES" OR NO 
У we eee eee ee ee ee 
ENTER THE ROW NUMBER OF THE ELEMENT TO BE CHANGED. 
2 
ENTER THE COLUMN NUMBER OF THE ELEMENT TO BE CHANGED. 
2 
THE ELEMENT Q( 2, 2)- 
0 
THE OUTPUT MEASUREMENT COST MATRIX ("О"-МАТВТХ)... 
1.00000 0.00000 
0.00000 0.00000 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO". 
Е . 
ENTER THE CONTROL COST WEIGHTING MATRIX {"R"-MATRIX} 
DIMENSION = # CONTROLS {NC} X # CONTROLS {NC}. 
YOU MAY ENTER THE MATRIX OR HAVE THE PROGRAM GENERATE 
A DESIRED MATRIX BY SELECTING ONE OF THE FOLLOWING: 
1 - NULL MATRIX 
2 - IDENTITY MATRIX 
3 - IDENTITY MATRIX SCALED BY A DESIRED VALUE 
4 - ANY OTHER MATRIX 
SELECT ANY OPTION 
4 


THE ELEMENT R( 1, 1)- 
.00002 
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iu CONTROPWS0ST MATRT1 SSS s т. | Am 
0.00002 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
ТУРЕ YES U OR "UNO 


FLAG/PARAMETER SETTINGS FOR THIS RUN ARE AS FOLLOWS: 


SYSTEM was input aS STABLE 


THE SAMPLE TIME INTERVAL FOR THE SYSTEM IS = 0.10 
DEBER OF SYSTEM - 2 

NUMBER OF CONTROLS - 1 

NUMBER OF OBSERVATIONS - 2 


NUMBER OF PROCESS NOISE SOURCES - 0 
OPTREG: DISCRETE OPTIMAL FILTER ANALYSIS...ORACLES PROGRAMS 
THE SAMPLE TIME INTERVAL FOR THE SYSTEM 15 - 0.10 

SYSTEM WAS INPUT AS STABLE 


А МАТКІХ 2 ROWS 2 COLUMNS 
1. 0000000Е+00 oO [55727202 
0.0000000Е%00 0. 3128365E 0] 

B MATRIX 2 ROWS 1 COLUMNS 
be 7 20959E-03 
6.3082559E-02 

H MATRIX 2 ROWS 2 COLUMNS 
1.0000000E+00 0.0000000E+00 
0.0000000E+00 1.0000000E+00 

8 МАТКІХ 2 ROWS 2 COLUMNS 
1. 0060000Е+00 0.0000000Е-00 
0.0000000Е%00 0.0000000Е>00 

R MATRIX 1 ROWS 1 COLUMNS 
2.0000000Е-05 


RICTNW: SUBROUTINE TO SOLVE DISCRETE STEADY-STATE 
RICCATI EQUATION BY THE NEWTON ALGORITHM. 


MATRIX (H TRANSPOSE)QH 


HTOH MATRIX 2 ROWS 2 COLUMNS 
1. 0000000Е+00 0.0000000Е%00 
0.0000000Е%00 0.0000000Е%00 


ІМ SUM, THE SEQUENCE OF PARTIAL SUMS HAS EXCEEDED STAGE 
50 without convergence 


FINAL VALUES OF P AND F AFTER 6 ITERATIONS TO CONVERGE 


р МАТКІХ 2 ROWS 2 COLUMNS 
1.6770747E+00 133590805 6E-07 
5 ))53850E-02 .9991775F 03 

F MATRIX 1 ROWS 2 COLUMNS 


1-0377 T9E+02 1.266 098E:01 


RESIDUAL ERROR IN RICCATI EQUATION 


EROR, MATRIX 2 ROWS 2 COLUMNS 
.3822048Е-09 8.4713316Е-10 
8 .4713316Е-10 9.72106 VESI] 
EIGENVALUES OF P 
EVLP MATRIX 2 ROWS 1 COLUMNS 
3.295 1535EFE-03 
1.6787757E+00 
CLOSED-LOOP RESPONSE MATRIX A-BF 
A-BF MATRIX 2 ROWS 2 COLUMNS 
6.2525884E-01 3.7153110E 02 
-6.952875в8Е 00 = 107720258500 


FOR THE ORIGINAL SAMPLED-DATA PROBLEM: 


FHAT MATRIX 1 ROWS 2 COLUMNS 
1.1037719F 02 1.2666098E+01 
2 COLUMNS 
1.1037719E+02 1.2666098E+01 
: E CONTROL GAINS CALCULATED.......O0PTREG COMPISEJBEDEENEC 
DO YOU WISH TO GENERATE AN "OPTMATD" DATA FILE TO SAVE THE: 
(A) (SYSTEM), (B) (CONTROL), 
{H} (OBSERVABLES), (G) (MEASUREMENT- FEED FWD. DIST.), 
(GAMD) (NOISE), (F) (CONTROL GAIN), 


(K) (FILTER GAIN), (Q) (OUTPUT COSTI 
(R) (CONTROL COST), (V1) (PROCESS NOISE INTENSITY (PSD)), 
(V2) (MEASUREMENT NOISE INTENSITY (PSD)), 
(S) (DISCRETE CROSS PRODUCT WEIGHTING) 
MATRICES FOR REENTRY TO THE OPTSYSX OR ORACLS PROGRAMS? 
(Y OR N) 
ORDSAL OPTIONS: 
1 -- KBFIL : DISCRETE KALMAN-BUCY FILTER ANALYSIS 
2 -- OPTREG: DISCRETE OPTIMAL REGULATOR ANALYSIS 
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ЕТ: RETURN ТО ОКАСІ,5 OPTIONS 
4 -- EXIT: RETURN TO ORACLS EXEC 
SELECT AN OPTION: 1, 2, 3, OR 4. 
NOTE: Кр ЕУ ANY TIME DURING ANY PROGRAM TWO SUCCESSIVE 


ENTERS" TERMINATES THE PROGRAM, SENDING YOU 
DEUS CTOLTHE ORACDLSUEXEC. 


— — — — е мм mmn m m cm cm cm cm cum ә. но оо оо мәз -. е ло мм ә. ә: е: -.  — н Әр  — s s s е е е е а мм мм е оо кз кз кз oom 


MESSAGE SUMMARY: MESSAGE NUMBER - COUNT 
187 & 


х THE ORACLS EXEC ALLOWS THE FOLLOWING OPTIONS: *%%2%% 
ШШ(СТаз БО5ТРАУ: ------- DISCRETE SYSTEM ANALYSIS------- 
(A-D CONVERSION, OPT. REG., KALMAN-BUCY 

FILTER AND DISCRETE TRANSIENT ANALYSIS) 

2 TRANFUNC FORTRAN: -DISCRETE TRANSFER FUNCTION ANALYSIS- 
(SYSTEM EIGENVALUES, OPEN LOOP, NOISE 
AND COMPENSATOR T.F.s, MODAL MATRICES) 


3 OPTPLOT FORTRAN: (PLOTTING OF TRANSIENT RESPONSES FOR 
STATES AND CONTROLS ) 


4 OPTGRAPH FORTRAN: ( POLE-ZERO, ROOT-LOCUS ) 

pe OULPUT SELECTION -- TERMINAL OR "A" DISK FILE 

6 HELP -- PROGRAM AND DATA FILE RELATIONSHIP EXPLANATIONS 
E LT -- RETURN TO CMS 

“у SELECT AN OPTION, ANY OTHER INPUT RETURNS YOU TO MENU ** 


Rea. 29/1. 


record off 
END RECORDING OF TERMINAL SESSION 


g 
> 
о 
39 
SNO 
O 
NO 
Э 
tri 
H 
Pr 
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APPENDIX D 
ORTRAL EXAMPLE RUN 


=» ар аә o ә. eee eee eee i — — ж. 


"HELP' INFORMATION OPTIONS: 
1 -- OPTSYSX-ORACLS; EQUATION & MATRICES NOTATION COMPARED 
2 -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
3 -- KBFIL: DISCRETE KALMAN-BUCY FILTER 
4 -- OPTREG: DISCRETE OPTIMAL REGULATOR 
5 -- ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
СШ ЕСТЕ RETURN TO ORACLS OPTIONS 
Е ооа: RETURN TO ORACLS EXEC 
SELEGPISN OPTION: 1,22 797 6. OR 7; 


еә сәр о. сәр ә. е е мм ме е е е се е е со е се се е е е е е се со е ее со ме оо ж-қ)» «ә» кк. «ко «ы» «кр «ы» «ар «кр ж.» «кы «ыы «кы «кр «к» ee Р 


ORTRAL IS A DISCRETE SYSTEM TRANSIENT аще РКОСКАМ 
ТНАТ WILL ACCEPT INPUT FROM ЕІТНЕК ТНЕ 5СКЕЕМ ОК ЕКОМ АМ 
ӘРІМАТО ат б IT USES THE DISCRETE MATRICES: {A}, 
ЕМ (H) ( G AND CALCULATES THE TRANSIENT RESPONSE 
Е А DISCRETE SYSTEM FOR A DESIRED NUMBER OF STAGES/POINTS. 
THE PROGRAM CALCULATES THE TRANSIENT RESPONSE OF A 15СКЕТЕ 
SYSTEM FOR A GIVEN SET OF INITIAL CONDITIONS; "Х(0)" 
BASED ON THE FOLLOWING EQUATIONS: 


X(I+1) = {A}*K(1) + (БУЗШ(Т) 
Y(I) - (H)*X(I) * (G)*U(I) 
UCI) - (ЕХ СБШ) 


NOTE: IF THE HATRIX eae T; ш STABLE THE 
STEADY STATE VALUE OF "X", VEN BY: 


X = ( (I) - (XA)-(B)(F)) ОЕ 


MAY BE COMPUTED AND PRINTED. ("-MATRIX INVERSE) 
NOTE: {UC} IS THE CONTROL INPUT (DRIVING FUNCTION) AND MAY 
BE SELECTED AS EITHER A "STEP" OR "RAMP'. 


ORTRAL OUTPUTS TO THE OPTPLOT DATA FILE THE "TIME", CONTROL 
VECTOR; "U", AND THE STATE VECTOR; "X" FOR PLOTTING VIA 
OPTPLOT FORTRAN AVAILABLE IN THE ORACLS EXEC. 


=» ар е ао е ес е е е см е се ес е е е се е е со мм е bo «нш» «ык» «кыр е ее со е е е е се со ое мо е ме е е е е е е е е е е е е е е е 


"HELP" INFORMATION OPTIONS: 
1 -- OPTSYSX-ORACLS; EQUATION & MATRICES NOTATION COMPARED 


2/2 


2 -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
3 -- KBFIL: DISCRETE KALMAN-BUCY FILTER 
4 -- OPTREG: DISCRETE OPTIMAL REGULATOR 
5 -- ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
fae] EXIT: RETURN TO ORACLS OPTIONS 
ЕСІ EXIT: RETURN TO ORACLS EXEC 
SELECT AN OPTION: 1,2, 3, &, 5, 6, OR 7 
6 
GENERAL ORACLS OPTIONS: 48  — 
l -- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
2. =] JORDSAL: “DISCRETE SYSTEM ANALYSIS 
3 -- ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
LN SHIELD: PROGRAM & EQUATION DESCRIPTIONS 
5 -- EX DE RETURN TO ORACLS EXEC 
SELECT AN OPTION: 1,2, 3, 4, OR 5. 
NOTE: SELECT OPTION - 5; IF YOU ARE USING THE OUTPUT TO 
DISK OPTION AND MOULD LIKE TO CHANGE FILENAMES 
NOTE: NORMALLY ANY TIME DURING THE PROGRAM TWO SUCCESSIVE 
NULL “ENTERS” TERMINATES THE PROGRAM, SENDING YOU 
BACK TO THE ORACLS EXEC 
3 
H МА" В tt AMD" "Н", " " 17 " 11 11 " " "R "T. ү?" AND 
"S" MATRICES FROM YOUR PREVIOUS LC MS WERE SAVED. 
THE FOLLOWING OPTIONS ARE AVAILABLE: 
l. USE ALL OF THE SAME MATRICES AGAIN. 
2. USE SELECTED MATRICES AGAIN. 
3. INPUT ALL NEW MATRICES. 
ЕМШЕК 1,2, ов 3. 
NOTE: EACH SAVED MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 
NOTE: A CHANGE IN THE ORDER OF THE SYSTEM; -NS- 
REQUIRES OPTION 3. 
2 


DO YOU WISH TO SAVE THE "A''-MATRIX FROM THE LAST 
RUN TO BE USED IN THE FOLLOWING RUN? 


NOTE: THE MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 

AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 
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TYPE "YES" OR "NO". 


У аас ЕСИ ИТИНИН 
DO YOU WISH TO SAVE THE "B'-MATRIX FROM THE LAST 
RUN TO BE USED IN THE FOLLOWING RUN? 
NOTE: THE MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 
TYPE "YES" OR "NO". 
ee es oo eee 
DO YOU WISH TO SAVE THE "'F"-MATRIX FROM THE LAST 
RUN TO BE USED IN THE FOLLOWING RUN? 
NOTE: THE MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 
TYPE "YES" OR "МО". 
п 
ENTER THE f OF PROCESS NOISE SOURCES (NG) OF THE 
( GAMD -MATRIX). 
THE DISCRETE SYSTEM EQUATION: 
о 
X(I+1) = {A}*X(I) + {B}*U(I) + {GAMD}*(WD+WO) 
{NS*NS ) {NS*NC } (NS*NG) 
NG ето 
0 
ЕМТЕК THE OF MEASUREMENTS OR OBSERVATIONS {NO} OF THE 
{"H"-MATRIX}. 
THE OUTPUT EQUATION: THE MEASUREMENT EQUATION: 
(NO) (NO) 
Y(I) = (H)*X(I)*(G)*U(I) 211) - (H)*X(I) + (G)*U(I) « V(I) 
(NO"NS)  (NO*NC) (NO*NS) (NO*NC) (NO*1) 
МС” 
0 


IS THE SYSTEM YOU WISH TO EVALUATE À STABLE SYSTEM ? 
(Y)ES OR (N)O ? 

NOTE: TO GET THE STEADY STATE VALUE IN "ORTRAL' 
YOU MUST ANSWER YES! 
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РЕ и КК жо жо -. -ә » ы о ә. әш 
= -- æ е е а е е е е е см е е е ао е е е е е е е е е е е е = «кы чик чик əс е -- а le 


THE SYSTEM MATRIX (A 'SMATRIK) D 
1.00000 0.08016 
00000 0.63128 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 


ТҮРЕ YES FOR NOME 


n 
THE CONTROL DISTRIBUTION MATRIX ("B'"-MATRIX)... 

198 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
Е TYPEM YES OR КОИ 

ENTER THE FEEDBACK GAIN CONTROL MATRIX {"F"-MATRIX}. 

DIMENSION = # CONTROLS (NC) X # STATES (NS). 
Мы ЕН Бе Тыр ES, TIS, "ORAS GENERATE 
1 - NULL MATRIX 
2 - ANY OTHER MATRIX 
SELECT AN OPTION 

2 
м, 
525.5 

THE ELEMENT F( 1, 2)- 
Шу. 23 
EE THE FEEDBACK GAIN CONTROL MATRIX ("F"-MATRIX) 

158.50000 17.33000 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 

т IPE "YES" OR КО“. 


мо ме мо ми е е ә» е мо чо ë e ә мо е | «шә е е е M ë е е е е е см е е е се е е е е ес е е е ее е е е е мл оо е мо е е е мо мо мо мо е е ми 


ENTER INITIAL CONDITIONS OF THE STATE: STATE VECTOR ("X(0)") 
DIMENSION - 3 OF STATES (NS) X T. 
SELECT THE APPROPRIATE INITIAL CONDITIONS: 
E ALL- I.C. S ARE ZERO 
2 - NON-ZERO I.C."S 
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ENTER: I OP 2 


2 
^. THE ELEMENT X( l, 1)= АШИ S s | 
z^ 
THE ELEMENT X( 2, 1)- 
0 
THE INITIAL CONDITIONS VECTOR ("X(0)").... 
0100000 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO". 
n 
HOW MANY POINTS DO YOU WANT THE SYSTEM EVALUATED AT? 
NETS =? (ENTER AN INTEGER ANSWER) 
3 


THE CONTROL LAW FOR ORACLS IS: 
U(I) = (F) X C T) t Kana) 
WHERE (UC) IS THE CONTROL INPUT (I.E. DRIVING FUNCTION) 
THE FOLLOWING DRIVING FUNCTIONS ARE AVAILABLE: 
1 - NO INPUT (I.E. (UC) IS A NULL VECTOR) 
2 - STEP INPUT 
3 - RAMP INPUT 
ENTER 0, 1, OR 2 FOR THE DESIRED DRIVING FUNCTION. 


2 

АТ WHAT STAGE DO YOU DESIRE THE STEP INPUT TO START? 
T(1) = ? 
0 


AT WHAT STAGE DO YOU DESIRE THE STEP INPUT TO STOP? 

T(2) = ° 

NOTE: TO GET THE STEADY STATE VALUE OF X(T) THE STOP STA 
MU INA (Т 


ST BE EQUAL OR GREATER THAN THE 
THE NUMBER OF POINTS.) 


GE 
L STAGE. n. 


-.-ә -. ме е мом е ме е е е е е е е е е е е е е е е е е е е е е е е е е е се се е е е е е е е е е е е е см см мо мо е е со мо е мо 


WHAT IS THE MAXIMUM VALUE OF THE STEP INPUT? 
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1 
DO YOU WISH TO DRIVE CONTROL # 1 ? 
ANSWER (Y)ES OR (N)O 
y 
FLAG/PARAMETER SETTINGS FOR THIS RUN ARE AS FOLLOWS: 
THE SAMPLE TIME INTERVAL FOR THE SYSTEM IS - 0.10000 
THE DRIVING FUNCTION IS A "STEP" INPUT WITH: 
START TIME - 0 
END TIME = 1 
МАХ. АМР. - 1.00000 


CONTROLS BEING DRIVEN: 
CONTROL # 1 -- YES 


ORDER OF SYSTEM = 2 


NUMBER OF CONTROLS = 1 
NUMBER OF OBSERVATIONS = 0 
NUMBER OF PROCESS NOISE SOURCES = 0 


-— ә = = еа е оа е оа е е е i i ee «мн | «ымы оа е «ж, 


THE SAMPLE TIME INTERVAL (DELT) FOR THIS SYSTEM IS = 0.1000 

THE DRIVING FUNCTION IS A "STEP" INPUT WITH: 
START TIME 
ND 


Е TIME 
MAX. AMP. 


CONTROLS BEING DRIVEN: 
CONTROL # 1 -- YES 


= = = = = нн «к. шн «шн. ат ее е е е е е е еш» оа аа е ат еа еа е е е е е е е оа е е «ә. «ма «ыма «ма еа е ое ое е е е ос е со осно е се е е е 


А МАТКІХ 2 ROWS 2 COLUMNS 
1.0000000E+00 9.0155729ËÉ 02 
0.0000000E+00 6.3129365E 01 
B MATRIX 2 ROWS 1 COLUMNS 
5.3950959E 05 
6.3082559Е-02 
H IS À NULL MATRIX 
G IS A NULL MATRIX 
2 COLUMNS 


F MATRIX 1 ROWS 
1.5850000Е+02 1.7330000E+01 
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STRUCTURE OF PRINTING TO FOLLOW 


STAGE TIME (SEC) = STAGE*DELT 
INPUT - UÇ(I) RANSPOSE; FOR APPROPRIATE DRIVING FUNCTION 
STATE - i ) TRANSPOSE; FROM X(I*1) - (A]*X(I) * (B)*U(I) 
CONTROL - U(I) TRANSPOSE; FROM U(I) “£ -(F)*x(1) * UC(I) 
0 0.00000 
UC(I): 0.1000000Е+01 
Х(І): 0.1000000Е%00 0.0000000Е+00 
U(I): -0.1585000E+02 
1 0.10000 
UC(I): 0.1000000Е+01 
Х(І): 0.4618773Е-01 -0.9998586Е»00 
U(I): 0.1100679Е+02 
2 0.20000 
UC(I): 0.0000000Е+00 
Х(І): 0.3412459Е-02 0.6314235Е-01 
0(І): -0.1635132E+01 
а 0.30000 
UC(I): 0.0000000E-*00 
X(I): 0.2922251Е-02 -0.6328755Е-01 
U(I): 0.6335965Е»00 
ко OPTPLOT DATE FILE GENERATED.....ORTRAL COMPLETE....... 
DO YOU WISH TO GENERATE AN "OPTMATD' DATA FILE TO SAVE THE: 
(A) (SYSTEM), (B) (CONTROL), 
(H) (OBSERVABLES), (G) (MEASUREMENT- FEED FWD. DIST.), 
(GAMD) (NOISE), (F) (CONTROL GAIN), 


(K) (FILTER GAIN), WW Q) (OUTPUT COST), 

(R) (CONTROL COST), (V1) (PROCESS NOISE INTENSITY (PSD)), 
(V2) (MEASUREMENT NOISE INTENSITY (PSD)), 

(S) (DISCRETE CROSS PRODUCT WEIGHTING) 

MATRICES FOR REENTRY TO THE OPTSYSX OR ORACLS PROGRAMS? 


GENERAL ORACLS OPTIONS: 
-- ORCONV: ANALOG TO DIGITAL MATRIX CONVERSION 
-- ORDSAL: DISCRETE SYSTEM ANALYSIS 
ORTRAL: TRANSIENT ANALYSIS OF A DISCRETE SYSTEM 
-- HELP: PROGRAM & EQUATION DESCRIPTIONS 
-- EXIT: | RETURN TO ORACLS EXEC 
SELECT AN OPTION: 1,2, 3, 4, OR 5. 


л г Wh P> 
1 
I 
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NOTE: SELECT OPTION - 5; IF YOU ARE USING THE OUTPUT TO 
DISK OPTION AND WOULD LIKE TO CHANGE FILENAMES 


NOTE: NORMALLY ANY TIME DURING THE PROGRAM TWO SUCCESSIVE 
NULL ENTERS' TERMINATES THE PROGRAM, SENDING YOU 
BACK TO THE ORACLS EXEC. 


-—— ас оо мис мо е е ее он мн е е е е eee е е он ме е е ем еее ы» 
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APPENDIX E 
TRANFUNC EXAMPLE RUN {1 (EIGENSYSTEM) 


BEGIN RECORDING OF TERMINAL SESSION 
R; T=0O.01/0.02 20:37 D 

controls 

DASD E LINKED R/O; R/W BY 0637P 

C (200) R/O 


PLOTTING IS DONE THRU DISSPLA. UP TO 204 OF YOUR DISK 
SPACE MAY BE NEEDED OTHERWISE ERROR MESSAGES 


m o m m + ata m^ m "to ale ale m m ale awe ale АА m < af. <t m ate ate m < < alo ale ale ale ala <t ale КА m m a! - m alo m - "o ata <t < m m m^ ale m m m + ata m m m ale ate 
0% 49% 6% 4% 4% 9% 47% S 9% 7% 4% 99 9% 9% 9% 9% 9% 94% 7% 9% 9% 7% 69% 4% 7% 7% 7% 06% 4% 4% ДО 4% 6% 79% 7% 49 7% 9% 9% 9% 9% 9% 9% 7% 7% 48S 7% 7% 9 9% 4% 49% 9 По 9% 9% 9% 7% 29 


CONTROLS CO* sISTS OF THREE INTERACTIVE PROGRAMS 


1 DACSAP SISO 
OPTSYS MIMO 
ORACLS MIMO 


HELP - PROGRAM DESCRIPTIONS 
SAMPLE PROBLEMS WITH PROGRAM DIRECTIONS 
ы DOCUMENTATION SOURCES 


ENTER^*T , 2 , SEA MO NUI 


Tat о наб аби аби аби а 0л абла e 


ale ale als alee ale alec ale a! оа ва ма steve we cet o ale steve we ve steve Seve 2.2.21... e se sese ses" uS im С este sles 4 Фа ФЬ ооа? e - ав em ta ata мама ма ма ња 
w 27% 07 9% 9% Zx 9% ом 07% 04% Фъ б 0% 9% 9% 4% 9% 9% FB 7% 7% 9 ge бо ZA 0*9 ж” 7% 7% оъ 9% 9% 9% 4% 


NOU UN 


NORMALLY IN ANY FORTRAN PROGRAM TWO NULL ENTREES TERMINATE 
THE PROGRAM SENDING YOU BACK TO THE CONTROLLING EXEC 


'251' REPLACES ' © (251) ° 2D ШШ "Ж 
F (231) R/O 
2517 REPLACES ' F (251) ' 
F (251) R/O 
хл ТНЕ ORACLS EXEC ALLOWS THE FOLLOWING OPTIONS: "ex 
l ORACLSK FORTRAN: — S DISCRETE SYSTEM. ANALYSIS E 
(A-D CONVERSION, OPT. RE KALMAN - BUCY 
FILTER AND DISCRETE TRANSIENT ANALYSIS) 
2 TRANFUNC FORTRAN: -DISCRETE TRANSFER FUNCTION ANALYSIS- 
(SYSTEM EIGENVALUES, OPEN LOOP, NOISE 
AND COMPENSATOR T.F:s, МОРАТ MATRICES} 


3. OPTPLOT FORTRAN: (PLOTTING OF TRANSIENT RESPONSES FOR 
STATES AND CONTROLS ) 


OPTGRAPH FORTRAN: ( POLE-ZERO, ROOT-LOCUS ) 

OUTPUT SELECTION -- TERMINAL OR “A” DISK FILE 

HELP -- PROGRAM AND DATA FILE RELATIONSHIP EXPLANATIONS 
ЕХІТ == ЕТО ОЗ 


< GN (д + 
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a SELECT AN OPTION, ANY OTHER INPUT RETURNS YOU TO MENU ** 


-. -. -ә = = 


“| о MW г (2 гә ғ 
I 


= = = -nmn мо е мн е сы әш сң ео әз å ә» -о ы» -ә өз ә» е» өза өң о ә ә» е» -. -. е е .. ә. ә. — i — — — «ж — ел «лә ел Ú «жә «ша «сө 


TRANFUNC OPTIONS: 
OPEN-LOOP EIGENSYSTEM ANALYSIS ONLY, (A) MATRIX 
OPEN-LOOP TRANSFER FUNCTION & EIGENSYSTEM 
NOISE TRANSFER FUNCTION & CLOSED-LOOP EIGENSYSTEM 
COMPENSATOR TRANSFER FUNCTION & ESTIMATOR EIGENSYSTEM 
MARKOV PARAMETER CHANGE 
HELP -- PROGRAM DESCRIPTION 
EXIT -- RETURN TO ORACLS ЕХЕС 


SELECT OPTION: 1, 2, 3, 4, 5, 6, OR 7 (SETS ITFOL FLAG) 


"HELP" PROVIDES AN OVERVIEW OF THE PROGRAM р 


NOTE: MODAL ANALYSIS AVAILABLE IN OPTIONS: 2, 3, 4 


= =ә = E Әәә а а е е е е е е е е мо =з тыз ==» тыз а е е е «ә». шз «ыш» «м». чө» - «м. ча» шз == «мә | «ә» «мн тыр ж =» «шн «м. | еш. ә аһ «иш» «мж» еш» «иш. == еш» «шә | еш» 
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AVAILABLE IN "TRANFUNC'" AND SOME OF THE EQUATIO 
USED IN THE CALCULATIONS. 


І -- TERMINAL VIEWING OF SELECTED INFORMATION 
2 -- TERMINAL VIEWING PLUS HARDCOPY TO YOUR "A" DISK 
3 -- EXIT: RETURN TO THE TRANFUNG OPTIONS 


NOTE: оон Pe. WILL PLACE THE FILE; “ТЕАМЕОМС НЕГР" 


YOUR "A" DISK CONTAINING THÉ INFORMATION YOU 
SELECTED TO VIEW AT THE TERMINAL. 


SELECT OPTION: TI 2, ОКУЗ 


=» m ee тыз тыз ae «еш» ек» ек» кыз жыз е е е е е е se se е см е е е е е е мо чє тз ж е е е е е е «а. ЕБ ЕБР ЕБР е е е е е е е е е со «шә «шы 


"HELP" INFORMATION OPTIONS: 


l -- TRANFUNC: EQUATIONS AND MATRICES DEFINED 

2 -- TRANFUNC: DISCRETE TRANSFER FUNCTION ANALYSIS 
Be = MODAL: DISCUSSION OF MODAL MATRICES 

4 -- EXIT: RETURN TO TRANFUNC OPTIONS 

БИ ЕН: RETURN TO ORACLS EXEC 


SELECT AN OPTION: 1, 2, 3, OR 4. 


= сәр ар ар аә сәр == ае е е == = е еш» Чыр е е е е е е е е е am am =» «шы. «к. «ы. чи» «к». е «ык. «к. «м» «к. е е мо е е е е се тыз == =» ©з е ое ое ë ‘M 


TRANFUNC IS A COMPLETELY INTERACTIVE DISCRETE TRANSFER 
FUNCTION ANALYSIS PROGRAM BASED ON THE FOLLOWING CONTROL 
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EQUATIONS: 
DISCRETE SYSTEM: (ORACLS NOTATION) 
SYSTEM EQUATION: 
X((I+1)%T) = (A)%K(I*T) + (B)%U(I*T) + (GAMD)*WD(I*T) 


OUTPUT EQUATION: Y(I) - (H)*X(I) * (G)*U(I) 
MEASUREMENT EQUATION:  Z(I) - (H)*X(I) * (G)*U(I) *V(I) 
CONTROL LAW: ШЕ- e 


COST EQUATION: aS GE 295 J: +) 


TCU (59 zutis TOS ati T5 Tor Т) 


WHERE ,T; IS THE SAMPLE TIME INTERVAL. 
WHER +" = TRANSPOSE. 
WHERE "g" IS LIMIT AS "N" APPROACHES INFINITY. 


DO YOU WISH A LIST OF MATRICES USED IN TRANFUNC? 
ANSWER (Y)ES OR (N) 


-. ә. --. «ж» «к» «к. «шә «к» «ж» «ше «кш» «ше «шә» “4 е см см е е е е е е е е е е е е е е ем е е се е е се е е е е е е е см «кш» «к» «к» ЧЧ е е е е «ке еш» “шы: | еш» 


CONTROL MATRIX: 

NOISE MATRIX: 

OBSERVABLES MATRIX: 

MEASUREMENT - FEED FWD. DIST. MATRIX: 
CONTROL GAIN MATRIX: 

FILTER GAIN MATRIX: 

OUTPUT COST MATRIX: 

CONTROL COST MATRIX: 

PROCESS NOISE MATRIX: 

MEASUREMENT NOISE MATRIX: 

CROSS PRODUCT WEIGHTING MATRIX 
RICCATI EQUATION SOLUTION MATRİX: 
NOISE VECTOR: 

CONSTANT NOISE VECTOR: 


SYSTEM MATRIX: tat 


ron 
5 
J 
I" 


куд о› << ДО АСО 


Ono КО 


"HELP" INFORMATION OPTIONS: 
-- TRANFUNC: EQUATIONS AND MATRICES DEFINED 
-- TRANFUNC: DISCRETE TRANSFER FUNCTION ANALYSIS 
MODAL: DISCUSSION OF MODAL MATRICES 
ccc RICE RETURN TO TRANFUNC OPTIONS 
== ШО RETURN TO ORACLS EXEC 


л + Ç N rr 
! 
I 


SELECT AN OPTION: 1, 2, 3, OR 4. 


TRANFUNC: DISCRETE TRANSFER FUNCTION ANALYSIS. 


102 


COMPUTES THE POLES, ZEROS, AND RESIDUES FOR A DESIRED 
TRANSFER FUNCTION. COMPUTATION OF MODAL MATRICES IS ALSO 
AVAILABLE DURING TRANSFER FUNCTION ANALYSIS. OUTPUTS 
RESULTS TO THE CORRESPONDING DATA FILES FOR PLOTTING VIA 
OPTGRAPH FORTRAN AVAILABLE IN THE ORACLS EXEC. OPTION 1 
MAY BE USED TO FIND THE EIGENVALUES OF THE SYSTEM MATRIX. 


DO YOU WISH TO VIEW THE TRANSFER FUNCTION EQUATIONS? 
ANSWER (Y)ES OR (N)O. 
TRANSFER FUNCTION EQUATIONS USED IN TRANFUNC 
--- OPEN-LOOP TRANSFER FUNCTION EQUATION: --- 
Z2/U * (H)*( z(I)-(A) )INV*(B) 
--- NOISE TRANSFER FUNCTION EQUATION: --- 
Z/U s (H)*( z(I)-(A)* (B) *(F) JINV*(GAMD) 
--- COMPENSATOR TRANSFER FUNCTION EQUATION: --- 
28 oe 2} Ae (By ee IK)%(H Y X ENV*TK) 


NOTE: THE COMPENSATOR TRANSFER FUNCTION IS FROM 
MEASUREMENT TO INPUT. 


-. - e с к á- á- оо мон е е ос е се е е е е е се зы. «жы мо е се е е е се е е ее е еле ее е е е е е е е е е е е ее еле е сы. = 


"HELP'' INFORMATION OPTIONS: 
1 -- TRANFUNC: EQUATIONS AND MATRICES DEFINED 


2 -- TRANFUNC: DISCRETE TRANSFER FUNCTION ANALYSIS 
-- MODAL: DISCUSSION OF MODAL HMATRICES 

4 -- EXIT: RETURN TO TRANFUNC OPTIONS 

КИС ЕХ RETURN TO ORACLS ЕХЕС 


SELECT AN OPTION: 1, 2, 3, ОВ 4. 


MODAL MATRICES MAY BE COMPUTED IN ALL TRANSFER FUNCTION 
ANALYSIS OPTIONS. MODAL MATRICES CAN BE USED IN DETERMINING 
WHETHER A SYSTEM IS CONTROLLABLE AND OBSERVABLE. IF THE 
MODAL DISTRIBUTION MATRIX, (TI*B) HAS NO ZERO ELEMENTS, THE 
SYSTEM IS COMPLETELY CONTROLLABLE. IF THE MODAL MEASUREMENT 
SCALING MATRIX, {H(BAR)*T} HAS NO ZERO ELEMENTS, THE 


103 


SYSTEM IS COMPLETELY OBSERVABLE/RECONSTRUCTIBLE. WITH THESE 
MATRICES THE USER CAN SEE AT À GLANCE WHETHER THE SYSTEM IS 
CONTROLLABLE AND OBSERVABLE WITH RESPECT TO ANY STATE/INPUT. 


SN ON Ui d^ C ГӘ = 
I 


"HELP' INFORMATION OPTIONS: 


1 -- ТКАМҒУМС: EQUATIONS AND MATRICES DEFINED 

2 -- TRANFUNC: DISCRETE TRANSFER FUNCTION ANALYSTS 
ЗБОРА: DISCUSSION OF MODAL MATRICES 

SNL ЕЛШІ RETURN TO TRANFUNC OPTIONS 
ЭББИ: RETURN TO ORACLS EXEC 


SELECT AN OPTION: T; 2, 3T ORG 


TRANFUNC OPTIONS: 
OPEN-LOOP EIGENSYSTEM ANALYSIS ONLY, {A} MATRIX 
OPEN-LOOP TRANSFER FUNCTION & EIGENSYSTEM 
NOISE TRANSFER FUNCTION & CLOSED-LOOP EIGENSYSTEM 
COMPENSATOR TRANSFER FUNCTION & ESTIMATOR EIGENSYSTEM 
MARKOV PARAMETER CHANGE 
HELP -- PROGRAM DESCRIPTION 
EXIT -- RETURN TO ORACLS EXEC 


SELECT OPTION: 1, 2, 3, 4, 5, 6, OR 7 (SETS ITFOL FLAG) 


THE 
MATRICES FROM YOUR PREVIOUS ЕЕ 


NOTE: MODAL ANALYSIS AVAILABLE IN OPTIONS: 2, 3, 4 


B "в" "GAMD' tt "VERS ML MA o E " " "р " VIS XD AND ке 
RUN W isum SAVED. 


THE FOLLOWING OPTIONS ARE AVAILABLE: 
1. USE ALL OF THE SAME MATRICES AGAIN. 
2. USE SELECTED MATRICES AGAIN. 

3. INPUT ALL NEW MATRICES. 


ЕМГ 2 FORS? 


NOTE: EACH SAVED MATRIX WILL BE REDISPLAYED AT 


THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 


NOTE: A CHANGE IN THE ORDER OF THE SYSTEM; -NS- 


REQUIRES OPTION 3 
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THE SYSTEM MATRIX ("A" - MATRIX)... 


0000 0.08016 

0000 0.63128 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
ТҮРЕ YES TORNO. 


FLAG/PARAMETER SETTINGS FOR THIS RUN ARE AS FOLLOWS: 


ITFOL IM IE ТТЕ1 ІТЕ2 ІТЕЗ IFGEW 
I 0 6 0 0 0 0 


THE SAMPLE TIME INTERVAL FOR THE SYSTEM IS = 0.10000 
ORDER OF SYSTEM = 2 

NUMBER OF CONTROLS = 1 

NUMBER OF OBSERVATIONS = 2 

NUMBER OF PROCESS NOISE SOURCES = 2 
ШАКТЕрМО:....... OPENT LOQEP ELGENE VALUE SANALYSIS ©... co ssw. ss 


TIROL IM TE ТТЕ] ТТЕ2 "ТЇТЕЗЛ,®#ТЕСЕ\ 
1 0 6 0 0 0 0 


THE SAMPLE TIME INTERVAL FOR THE SYSTEM IS = 0.10000 
А MATRIX 2 ROWS 2 COLUMNS 
1.0000000Е%00 "01557 2 5E-02 
0.0000000E+00 009120990 9JE - Om 
IPPENSISUUPCESGENVALUEST-....... DET(ZI-A).. 


ШООООООЕ+ОО Т 6.3128&Е-01: 


ОРЕМ LOOP RIGHT EIGENVECTOR MATRIX....... DT 
1.000000D+00 -2.124296D-01 
0.000000D+00 9.771764D-01 
OPEN LOOP LEFT EIGENVECTOR MATRIX...... Т-ІМУ 
1.000000р+00 2.173913р-01 
0.000000р+00 1.023357р+00 
 TRANFUNG COMPLETED ......- FILES FOR OPTGRAPH GENERATED.. 
DO YOU WISH TO GENERATE AN "OPTMATD" DATA FILE TO SAVE THE: 
(A) (SYSTEM), (B) (CONTROL), 
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{Н} (OBSERVABLES), {G} (MEASUREMENT- FEED FWD. DIST.), 


(GAMD) (NOISE), (F) (CONTROL GAIN), 
(K) (FILTER GAIN), (Q) (OUTPUT COST), 
(R) (CONTROL COST), (V1) (PROCESS NOISE INTENSITY (PSD)) 


(V2) (MEASUREMENT NOISE INTENSITY (PSD)), 
(S) (DISCRETE CROSS PRODUCT WEIGHTING), 
MATRICES FOR REENTRY TO THE OPTSYSX OR ORACLS PROGRAMS? 
(Y OR N) 


хеее THE ORACLS EXEC ALLOWS THE FOLLOWING OPTIONS: #34 


1 ORACLSX FORTRAN: TEER E DISCRETE SYSTEM ANALYSIS------- 
(A-D CONVERSION, OPT. REG., KALMAN-BUCY 
FILTER AND DISCRETE TRANSIENT ANALYSIS) 


2 TRANFUNC FORTRAN: -DISCRETE TRANSFER FUNCTION ANALYSIS- 
(SYSTEM EIGENVALUES, OPEN LOOP, NOISE 
AND COMPENSATOR T.F.s, MODAL MÀTRICES) 


3  OPTPLOT FORTRAN: (PLOTTING OF TRANSIENT RESPONSES FOR 
STATES AND CONTROLS) 


4 OPTGRAPH FORTRAN: ( POLE-ZERO, ROOT-LOCUS ) 

5 OUTPUT SELECTION -- TERMINAL OR "A" DISK FILE 

6 HELP -- PROGRAM AND DATA FILE RELATIONSHIP EXPLANATIONS 
77 "EXIIT --“RETURN“TOTCHS 

** SELECT AN OPTION, ANY OTHER INPUT RETURNS YOU TO MENU ** 
7 


DASD 200_DETACHED 
R; D20.27/ 1929. 16:22:67 


record off 
END RECORDING OF TERMINAL SESSION 
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APPENDIX F 
TRANFUNC EXAMPLE RUN #2 (OPEN-LOOP TF) 


TRANFUNC OPTIONS: 


- OPEN-LOOP EIGENSYSTEM ANALYSIS ONLY, {A} MATRIX 

- OPEN-LOOP TRANSFER FUNCTION & EIGENSYSTEM 

- NOISE TRANSFER FUNCTION & CLOSED-LOOP EIGENSYSTEM 
COMPENSATOR TRANSFER FUNCTION & ESTIMATOR EIGENSYSTEM 
- MARKOV PARAMETER CHANGE 

- HELP -- PROGRAM DESCRIPTION 

- МЕНИ сс RETURN TO ORACLS EXEC 


EEEECT.OPTTON: І, 2, 3, 4, 5, 6, OR 7 (SETS ITFOL FLAG) 
NOTE: M@DAL ANALYSIS AVAILABLE IN OPTIONS: 2, 3, 4 


OPEN-LOOP TRANSFER FUNCTION OPTIONS: 
L OREHI) {A} YINV=(B) 


OPTTION L -- PORES, RESIDUES, AND ZEROS COMPUTED 
OPTION 2 -- ONLY POLES AND ZEROS COMPUTED 
OPTION 3 -- ONLY POLES AND RESIDUES COMPUTED 
OPTION 4 -- EXIT -- RETURN TO TRANFUNC OPTIONS 


SELECT OPTION: 1, 2, 3, ОК 4 fonts = 1,2 ,06,3) 
NOTE: POLES AND ZEROS MUST BE SELECTED IF YOU 


TO USE OPTGRAPH 


NOTE: REQUIRES МС>0, МО>0 


1 
ТНЕ MALE "з САМР tt "H' Ec и ! R" V " MOM AND MEG TT 
MATRICES Mog YOUR PREVIOUS RE NOE d SAVED. 


NOPE: 


THE FOLLOWING OPTIONS ARE AVAILABLE: 

1. USE ALL OF THE SAME MATRICES AGAIN. 
2. USE SELECTED MATRICES AGAIN. 

3. INPUT ALL NEW MATRICES. 

Peek Ly 2,20R =o". 

EACH SAVED MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 
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AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 


NOTE: A CHANGE IN THE ORDER OF THE SYSTEM; -NS- 
REQUIRES OPTION 3. 


ғ. о. іш» om oem ee е е сәз UO әәә Әәә әәә мо . ме ее е е ее е мо мн «ә. е е е «м» е ее е m е см е есен е «лы тə» т» ты» зә» «зә «ә» е тə» сә» ел o o — umbo um 


THE SYSTEM MATRIX {'A'-MATRIX}... 


000 0.08016 
000 0.63128 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
ТҮРЕ "YES" OR NO. 
THE CONTROL DISTRIBUTION MATRIX ("B'-MATRIX)... 


0.00340 
0.06308 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO". 
THE MEASUREMENT SCALING MATRIX ("H'-MATRIX)... 


1.00000 0.00000 
0.00000 1.00000 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO". 


n 
DO YOU WISH TO DETERMINE THE MODAL DISTRIBUTION 
AND GAIN MATRICES? 
TYPE (Y)ES OR (N)O 
у 


FLAG/PARAMETER SETTINGS FOR THIS RUN ARE AS FOLLOWS: 


ITFOL IM IE ITFl ITF2 ITF3 IFGFW 
2 l 6 1 0 0 0 


THE SAMPLE TIME INTERVAL FOR THE SYSTEM IS = 0.10000 
ORDER OF SYSTEM - 2 

NUMBER OF CONTROLS = 1 

NUMBER OF OBSERVATIONS = 2 

NUMBER OF PROCESS NOISE SOURCES = 2 

IRANRBÜUNG ec OPEN LOOP TRANSFER FUNCTION ANALYSIS....... 
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ШЕНІ ІМ ТЕ ТТЕ1 ІТЕ2 ІТЕЗ ТЕСЕМ 
2 1 6 r О 0 0 
THE SAMPLE TIME INTERVAL FOR THE SYSTEM IS = 0.100 
À MATRIX 2 ROWS 2 COLUMNS 
1.0000000Е+00 5.015572% ue 
0.0000000E+00 03126365 EN 
B MATRIX 2 ROWS 1 COLUMNS 
2250250959 E-03 
8893052 559E-02 
H MATRIX 2 ROWS 2 COLUMNS 
1.0000000Е%00 О. 0000000Е+00 
О. 0000000Е+00 1. 0000000Е+00 
ШЕБПИБФОРИЕТСЕХММАВОЕВ.....2.... DET(ZI-A).. 
1.00000E+00: 6.31284E-01: 
OPEN COOP RICHT EIGENVECTOR MATRIX....... Toa 
1.000000D+00 -2.124296D-01 
0.000000D+00 9.771764D-01 
OFPENSEOOPMEEFT EIGENVECTOR MATRIX TI. T= INV. 
1.000000D+00 2.173913D-01 
0.000000D+00 1.023357D+00 
ИИ а а МОЮАТ, ЖИИ Т@Е5. о оо. а 
MODAL CONTROL DISTRIBUTION MATRIX...... В 


lee LOS 
6.4555 


OD-02 


7 
96D-02 


MODAL MEASUREMENT SCALING MATRIX...H(BAR)*T.. 


JOD - 0! 
64D-01 


0р+00 
OD+00 


1.000 


00 -2.1242 
0.00000 17 


du ud 


TF FOR INPUT NO. 1 AND OUTPUT NO. 
ORDER OF NUMERATOR 1 
NUMERATOR EIGENVALUES (INCLUDING 

( 0.000000)+J( 0.000000) 
(-0.858049)«J( 0.000000) 


RESTDUESTAT THE POLES: 
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TF GAIN 


TL: 


0.3395E-02 
EXTRANEOUS ZERO VALUES): 


P OL ES 
REAL(A) IMAG(B) 


( 1.000000)«J( 0.000000) 
( 0.631284)«J( 0.000000) 


КЕ 5 I D USERS 


( 0.017109)*EXP(A*T) 
(-0.013714)*EXP(A*T) 


TF FOR INPUT NO. l1 AND OUTPUT NO. 2: 


ORDER OF NUMERATOR - 1 


TF GAIN = 0.6308Е-01 


NUMERATOR EIGENVALUES (INCLUDING EXTRANEOUS ZERO VALUES): 


( 1.000000)+J( 0.000000) 
( 0.000000)+J( 0.000000) 


RESIDUES AT THE POLES: 


P O DE S 
REAL(A) IMAG(B) 


( 1.000000)+J( 0.000000) 
( 0.631284)+J( 0.000000) 


— — = == = == ек» am em «ж. ча = ew ow = ow = = owe = = еш» = aw = ae шыт 


.IRANEUNG COMPLETED. ..... 


DO YOU WISH TO GENERATE AN 


КЕ 5 МЕ S 


( 0.000000)*EXP(A*T) 

( 0.063083)*EXP(A*T) 
..FILES FOR OPTGRAPH GENERATED. . 
"OPTMATD" DATA FILE TO SAVE THE: 


(A) (SYSTEM), (B) (CONTROL), 

(H) (OBSERVABLES), {6} (MEASUREMENT- FEED FWD. DIST.), 
(GAMD) (NOISE), (F) (CONTROL GAIN), 

(K) (FILTER GAIN), {Q} (OUTPUT COST), 

{R} (CONTROL COST), (Vl) (PROCESS NOISE INTENSITY (PSD)) 


(V2) (MEASUREMENT NOISE INTENSITY (PSD)), 
(5% (DISCRETE CROSS PRODUCT WEIGHTING), 
MATRICES FOR REENTRY TO THE OPTSYSX OR ORACLS PROGRAMS? 
(Y OR N) 
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APPENDIX G 
TRANFUNC EXAMPLE RUN #3 (NOISE TF) 


= мо мо мн е á- е мо е ес мо мм мо ӘӘ ә Әә Әлә е е ее е е е е ое е се мо мо е ем а мм е е е ом ме еш» еш» өше  — жә өш» 


TRANFUNC OPTIONS: 
- OPEN-LOOP EIGENSYSTEM ANALYSIS ONLY, (A) MATRIX 
- OPEN-LOOP TRANSFER FUNCTION & EIGENSYSTEM 
- NOISE TRANSFER FUNCTION & CLOSED-LOOP EIGENSYSTEM 
COMPENSATOR TRANSFER FUNCTION & ESTIMATOR EIGENSYSTEM 
- MARKOV PARAMETER CHANGE 
- HELP -- PROGRAM DESCRIPTION 
- EXIT -- RETURN TO ORACLS EXEC 
ERECT OPTION: 1 2,3 4, 5 6. ORVZ. (SETS ITFOL FLAG) 
NOTE: MODAL ANALYSIS AVAILABLE IN OPTIONS: 2, 3, 4 


ч ON Ui AUN El 
! 


= ею сә иә е е е е е е е е е е е е е е е е е е оз е е; оо е е е е е е е е е е е е е е е е е е е sw оо е ss е е ss е со е е е е == 


NOISE TRANSFER FUNCTION OPTIONS: 
ieee { 2{т}={А}+{В}%З{Е} JINV™{GAMD)} 
OPTION 1 -- POLES, RESIDUES, AND ZEROS COMPUTED 
OPTION 2 -- ONLY POLES AND ZEROS COMPUTED 
OPTION 3 -- ONLY POLES AND RESIDUES COMPUTED 
OPTION 4 -- EXIT -- RETURN TO TRANFUNC OPTIONS 

SELECT OPTION 1, 2, 3, OR 4 (SEES TLE? = 12,3) 
NOTE: NOISE TF FUNCTION THRU CLOSED LOOP SYSTEM; 
NOTE: POLES AND ZEROS MUST BE SELECTED IF YOU 


ish to use OPTGRAPH. 
NOTE: REQUIRES NC>0, NO>0, NG>O 


ма мо а е е ес е ee ee ee ee е е е е е е е оо е е оо е е со е е е е е е 


ТНЕ Шыда "B" "САМЮ:' tt tt e t р 1! “К 11 tt үл? д! AND EST 
MATRICES FROM YOUR PRÉVIOUS ORACLE iu S d SAVED. 


THE FOLLOWING OPTIONS ARE AVAILABLE: 
l. USE ALL OF THE SAME MATRICES AGAIN. 
2. USE SELECTED MATRICES AGAIN. 

3. INPUT ALL NEW MATRICES. 


ENTER IM 2, “ORS: 
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NOTE: EACH SAVED MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 


NOTE: A CHANGE IN THE ORDER OF THE SYSTEM; -NS- 
REQUIRES OPTION 3. 


1 
THE SYSTEM MATRIX (A MATRI 
1.00000 0.08016 
0.00000 0.63128 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE- YES OR NOL: 
n 
THE CONTROL DISTRIBUTION MATRIX {"B"-MATRIX}... 
0.00340 
0.06308 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE "YES" OR "NO". 
n 


THE PROCESS NOISE DISTRIBUTION MATRIX ("GAMD"-MATRIX)... 

0 0.00000 

0 1.00000 

DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE YES OR S NO 


n 
THE MEASUREMENT SCALING MATRIX {"H"-MATRIX}... 
180008 2:8888 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
Е TYPE "YES" OR "NO". 
THE FEEDBACK GAIN CONTROL MATRIX {"F"-MATRIX} 
110.40000 12.66000 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
И TYPE "YES" OR "NO". 


DO YOU WISH TO DETERMINE THE MODAL DISTRIBUTION 
AND GAIN MATRICES? 


TYPE (Y)ES OR (N)O 
FLAG/PARAMETER SETTINGS FOR THIS RUN ARE AS FOLLOWS: 
ІТЕО, ІМ ІЕ ІТЕІ ІТЕ2 ІТЕЗ IFGFW 
3 1 6 0 1 0 0 


THE SAMPLE TIME INTERVAL FOR THE SYSTEM IS = 0.10000 
ORDER OF SYSTEM = 2 

NUMBER OF CONTROLS - 1 

NUMBER OF OBSERVATIONS = 2 

NUMBER OF PROCESS NOISE SOURCES = 2 
БШЕАКЕШМС:............ NOISE TRANSFER FUNCTION ANALYSIS...... 


FTFOL ІМ TE ІТЕ! ІТЕ2 ІТЕЗ ТЕСЕУ 
3 il 6 0 1 0 0 


IHE SAMPLE TIME INTERVAL FOR IHE SYSTEM IS - 0.10000 
А МАТКІХ 2 ROWS 2 COLUMNS 
1.0000000E-*00 20155722 02 
0.0000000Е%00 5. 31283656501 
B MATRIX 2 ROWS 1 COLUMNS 
ВЕ 9309539E-93 
53082559 E02 
САМО MATRIX 2 ROWS 2 COLUMNS 
1. О000000Е+00 0.0000000E-*00 
О. 0000000Е+00 1. 0000000Е+00 
Н МАТКІХ 2 ROWS 2 COLUMNS 
1.0000000Е%00 0. 0000000Е+00 
0.0000000Е%00 1. 0000000Е+00 
F MATRIX 1 ROWS 2 COLUMNS 
1.1040000Е+02 1.2660000E+01 
OPENGL OOFMETGENVALUES: -ap t aeaa DET(ZI-A).. 
1.00000E+00: 6.31284E-01: 
OPEN LOOP RIGHT EIGENVECTOR MAIRIX....... Т 
1.000000D+00 -2.124296D-01 
0.000000D+00 9.771764D-01 
OPENTEOOP LEFT EIGENVECTOR MATRIX.. ..: ПЕЕ 
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MODAL CONTROL DISTRIBUTION MATRIX. ..... TB 


1.710870D-02 
6.455596D-02 


MODAL PROCESS NOISE DISTRIBUTION MATRIX...TI*GAMD.. 


1.000000D+00 2.173913D-01 

0.000000D+00 1.023357р»00 

MODAL MEASUREMENT SCALING MATRIX...H(BAR)*T.. 
1.000000р+00 -2.124296р-01 

0.0000000%00 9.771764D-01 

THE MODAL CONTROL GAINS.: SE PT.. 
1.104000D+02 -1.108118D+01 

.. IHE CLOSED LOOPSDYNAMICS MATRIX ASE p m 
6.251814D-O01, 3.717381D-02 
-6.964315D+00 -1.673415D-01 


CLOSED LOOP REGULATOR E-VALUES..DET(ZI-A*B*E).- 
2:29920E-0 ME 166FE-01: 


C-LOOP RIGHT EIGENVECTOR MADRE T ни а dos. 
-5.689885D-02 -4.582878D-02 

1.0000000%00 0.000000D-00 

CONTROL EIGENVECTOR MATRIX 0 PSS esa 
C-LOOP ОРТ. КЕС. LEFT E-VECTOR MATRIX M TNV 
0.000000D+00 1.000000D+00 
-2 .182035D+01 -1.241553 D 00 


.NOISE TRANSFER FUNCTIONS THROUGH THE CLOSED-LOOP SYSTEM.. 


TF FOR INPUT NO. 1 AND OUTPUT NO. 1: 
ORDER OF NUMERATOR = 1 TF GAIN = 0.1000E+01 
NUMERATOR EIGENVALUES (INCLUDING EXTRANEOUS ZERO VALUES): 
(-0.167342)+J( 0.000000) 
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( 0.000000)+J( 0.000000) 


RESIDUES AT THE POLES: 


ПОЕ РЕ ТЕ 5 
ВЕАҺ(А) ІМАС(В) 


( 0.228920)-J( 0.319166) ( 1.000000) EXP(A*T)*COS(B*T) 
( 0.228920)+J(-0.319166) ( 1.241553) EXP(A*T)*SIn(B*T) 


EEOEORSENPUT NO. TI AND OUTPUT NO. 2: 
NO FINITE ZEROS. TF GAIN = -0.6964E+01 


RESIDUES AT THE POLES: 


POLES PMES I DUES 
REAL(A) IMAG(B) 


( 0.228920)+J( 0.319166) ( 0.000000) EXP(A*T)*COS(B*T) 
( 0.228920)5Ј(-0.319166)  (-21.820348) EXP(A*T)*SIN(B*T) 


CET FOR INPUT NO. 2 AND OUTPUT NO. 1: 
NO FINITE ZEROS. ТЕ САҒЫ = 0.371/7Е-01 


FESIDUES AT THE POLES: 


POLES RESIDUES 
REAL (A) IMAG (B) 


( 0.228920)+J( 0.319166) ( 0.000000) EXP(A*T)*COS(B*T) 
( 0.228920)*J(-0.319166) ( 0.116472) EXP(A*T)*SIN*B*T) 


TF FOR INPUT NO. 2 AND OUTPUT NO. 2: 
ORDER OF NUMERATOR = 1 TF GAIN = 0.1000Е+01 
NUMERATOR EIGENVALUES (INCLUDING EXTRANEOUS ZERO VALUES): 
( 0.625181)+J( 0.000000) 
( 0.000000) +7( 0.000000) 


RESIDUES AT THE POLES: 


POL DELE I p E S 

REAL(A) IMAG (B) Mo | 
( 0.228920)>7( 0.319166) ( 1.000000) EXP(A*T)*COS(B*T) 
( 0.228920)+3(-0.319166) (-1.241553) EXP(A*T)*SIN(B*T) 


.- ТКАМЕШЫС COMPLETEDE a. FILES FOR OPTGRAPH GENERATED.. 
DO YOU WISH TO GENERATE AN "OPTMATD" DATA FILE TO SAVE THE: 


(A) (SYSTEM), (B) (CONTROL), 

(H) (OBSERVABLES), {6} (МЕАЅОВЕМЕМТ- FEED FWD. DIST.), 
{GAMD} (NOISE), {F} (CONTROL GAIN), 

(к) (FILTER GAIN), (Q) (OUTPUT COST), 

(В) (CONTROL COST), (V1) (PROCESS NOISE INTENSITY (PSD)) 


(V2) (MEASUREMENT NOISE INTENSITY (PSD)), 
(S) (DISCRETE CROSS PRODUCT WEIGHTING), 
MATRICES FOR REENTRY TO THE OPTSYSX OR ORACLS PROGRAMS? 
(Y OR N) 


-— = ес ес е а е сыз тыз тыз е сыз тыз жыр жыр сыз сыз тыз ЧЕР чыз е е е е е е е е е е е сыр е е е е е е е е е се е е е е е е е е е е е есен е е ол 


еә сәр әәә әәә е е е е е е е ее е е е е е е е ее е е е е е е е = е е е е ө е е се ее е = е е е = е е е е е е е е е е е е е 
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APPENDIX H 
TRANFUNC EXAMPLE RUN #& (COMPENSATOR TF) 


TRANFUNC OPTIONS: 
OPEN-LOOP EIGENSYSTEM ANALYSIS ONLY, (A) MATRIX 
OPEN-LOOP TRANSFER FUNCTION & EIGENSYSTEM 
NOISE TRANSFER FUNCTION & CLOSED-LOOP EIGENSYSTEM 
COMPENSATOR TRANSFER FUNCTION & ESTIMATOR EIGENSYSTEM 
MARKOV PARAMETER CHANGE 
HELP -- PROGRAM DESCRIPTION 
EXIT -- RETURN TO ORACLS EXEC 


ELECT OPTION. 1, 2, 3, 4, 5,96, OR-7 (SETS ITFOL FLAG) 


THE 


NOTE: MODAL ANALYSIS AVAILABLE IN OPTIONS: 2, 3, 4 


=» MEE ME NN ©} ое е ee е е е е е е е е е е мн мо ше мо әнә әнә мо мом 


COMPENSATOR TRANSFER FUNCTION OPTIONS: 


G О О (F) OE (T) Ar B) (F) K) (H) }INV=(K} 
OPTION 1 -- POLES, RESIDUES, AND ZEROS COMPUTED 
OPTION 2 -- ONLY POLES AND ZEROS COMPUTED 
OPTION 3 -- ONLY POLES AND RESIDUES COMPUTED 
OPTION 4 -- EXIT -- RETURN TO TRANFUNC OPTIONS 
SELECT OPTION 1, 2, 3, OR 4 (nsn 5 1.2.3) 


NOTE: COMPENSATOR T.F. IS FROM MEASUREMENT TO INPUT 


NOTE: A COMPENSATOR TRANSFER FUNCTION MAY BE COMPUTED 
ONLY IF BOTH THE REGULATOR ДЕ) AND FILTER (K) 
GAINS ARE AVAILABLE. THEY MAY BE CALCULATED IN 
THE "OPTREG" AND "KBFIL" OPTIONS OF ORACLSX. 

NOTE: POLES AND ZEROS MUST BE SELECTED IF YOU WISH 
TO USE OPTGRAPH. 

NOTE: REQUIRES NC»0, NO>O 


чат ме = ат ew eS е е е е е се см см мо е се се е со со ое со ее со со см ем осмо е мо е со см е со е смо е се е со со со мо е е -- 


=— = = = es = ee et o GRO RO ее е е он ено а а а ақы а е е е се е е се со е е е е е е он со см он е е е е е е е со со мо оо m 


EA "Б "GAMD' "Hg" п: Е DK " R" EVE MO AND С! t 


MATRICES FROM YOUR PRÉVIOUS Qus B WERE SAVED. 


THE FOLLOWING OPTIONS ARE AVAILABLE: 
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1. USE ALL OF THE SAME MATRICES AGAIN. 
2. USE SELECTED MATRIGES “AGAIN. 
3. INPUT ALL NEW MATRICES. 


ENTER 1, 2, OR 3. 


NOTE: EACH SAVED MATRIX WILL BE REDISPLAYED AT 
THE PROPER INPUT SEQUENCE INTERVAL 
AND YOU WILL HAVE THE OPTION OF CHANGING 
INDIVIDUAL MATRIX ELEMENTS. 


NOTE: A CHANGE IN THE ORDER OF THE SYSTEM S MS. 
REQUIRES OPTION 3. 


ІВ 
THE SYSTEM MATRIX ("A"-MATRIX)... 
0:00000 0:63128 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
З TYPE "MES" ORO 
THE CONTROL DISTRIBUTION MATRIX {"'B"-MATRIX}... 
0106308 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
x TYPE "NES “OR™ to = 
THE MEASUREMENT SCALING MATRIX ("H'-MATRIX)... 
0:00000 1:00000 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
M ТҮРЕ YES ОРЕ МО 
THE FEEDBACK GAIN CONTROL MATRIX ("F'"-MATRIX) 
110.40000 12.66000 
DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
Е ТҮРЕ YES "OR NO ` 


THE FEEDBACK GAIN ESTIMATOR MATRIX {"K"-MATRIX} 


0.61914 0.04938 
0.00338 0.34673 


DO YOU WISH TO CHANGE THE VALUE OF ANY MATRIX ELEMENT? 
TYPE YES? (OR NO ` 
Пп 
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DO YOU WISH TO DETERMINE THE HODAL DISTRIBUTION 
AND GAIN MATRICES? 


TYPE (Y)ES OR (N)O 


т“ ме ма че ма = е е се е мм е см е е е е е е е е е е е е 1l е е е к. lcT ll MI шо шо шо ош 


FLAG/PARAMETER SETTINGS FOR THIS RUN ARE AS FOLLOWS: 


ITFOL ІМ IE ITF1 ITF2 ІТЕЗ ТЕСЕМ 
4 1 6 0 0 1 0 


ШИЕ SAMPLE TIME INTERVAL FOR THE SYSTEM IS 
ORDER OF SYSTEM 2 

NUMBER OF CONTROLS 
NUMBER OF OBSERVATIONS 2 
NUMBER OF PROCESS NOISE SOURCES 


0.10000 


— 
— 


1 


IBESSNEONC:...... COMPENSATOR TRANSFER FUNCTION ANALYSIS...... 
ITFOL ІМ ТЕ ІТЕІ! ITF2 ITF3  IFGEW 
4 l 6 0 0 il 0 
THE SAMPLE TIME INTERVAL FOR THE SYSTEM IS = 0.10000 
А МАТКІХ 2 ROWS 2 COLUMNS 
1.0000000E+00 SUELOS 
0.0000000Е%00 CIUS OST OE 
B MATRIX 2 ROWS 1 COLUMNS 
252950959505 
03082559E-02 
H MATRIX 2 ROWS 2 COLUMNS 
1.0000000E+00 0.0000000E+00 
0.0000000E+00 1.0000000E+00 
F MATRIX 1 ROWS 2 COLUMNS 
1.1040000E+02 1.2660000E+01 
K MATRIX 2 ROWS 2 COLUMNS 
6. 1913667/E-01 4.9382360E-02 
ЭЛЕ ӘСЕ - 03 3.4673134E-01 
OPPP SPO OP EIGENVALUES.......... БЕШЕ Ај. 
1.00000Е%00: 6.31284Е-01: 
OPEN LOOP RIGHT EIGENVECTOR MATRIX....... ТЕ 
1.000000D+00 -2.124296D-Q1 
0.000000р-00 9.771764D-01 
OPEN E@OP LEFT ELGENVECTOR MATRIX...... T-INV 


1.000000D+00 


ЖЕТТІ ВИ 


9 


0.000000D+00 1.023357D+00 


2 
А 


MODAL CONTROL DISTRIBUTION MATRIX...... TI p: 

1.710870D-02 

6.455596D-02 

MODAL MEASUREMENT SCALING MATRIX...H(BAR)*T.. 
1.0000000-00 -2 .14242960- И 

0.0000000+00 29.771764D-01 

THESMODAL CONTROENGAINS....:2 21201227 2 E DOM. 


1.104000р+02 -1.108118р+01 
MODAL FILTER STEADY STATE GAINS........ TI*K.. 


6.198718D-01 ы 247587D-Q1 
3.460760D-03 .548298D-01 


CLOSED LOOP ESTIMATOR E-VALUES..DET(ZI-A*K*H).. 


.61982E+00: 9277535905 01: 


LOOP RIGHT EIGENVECTOR MATRIX. ........ Е 


986 1D-01 im ШЕ 01 
90930-03 9.802334р-01 


(т 

9 

6 

MEASUREMENT EIGENVECTOR MATRIX..... Н(ВАК УЕ 

9.9998610-01 -1.978447D-01 

5.269093D-03 9.80233383D-01 

-LOOP OPT. FILTER І.ЕЕТ Е-УЕСТОК ИР ie 

5D-01. 2.016226D 0: 

9D-03 1.019081D+00 

COMPENSATOR TRANSFER FUNCTIONS FROM MEAS. TO INPUT 
U/Z = -F*(zI-A+B*F+K*H)INV*K 
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THE CLOSED LOOP FILTER DYNAMICS MATRIX ...{A-KH}.. 


TF DENOMINATOR EIGENVALUES: 
(uty T1854 )+J( 0.564795) 
( 0.711854) +7(-0.564795) 


TF FOR INPUT NO. 1 AND OUTPUT NO. 1: 
ORDER OF NUMERATOR = 1 TF GAIN = 0.6840E+02 


NUMERATOR EIGENVALUES (INCLUDING EXTRANEOUS ZERO VALUES): 
( 0.977322)+J( 0.000000) 
( 0.000000)+J( 0.000000) 


RESIDUES AT THE POLES: 


POLL ES PE SIDUES 
REAL (A) ІМАС(В) 


( 0.711854)«J( 0.564795)  ( 68.395501) EXP(A*T)*COS(B*T) 
( 0.711854)*J(-0.564795)  (-32.147575) EXP(A*T)*SIN(B*T) 


TF FOR INPUT NO. 2 AND OUTPUT NO. 1: 
ORDER OF NUMERATOR - 1 TF GAIN = 0.9841E+01 
NUMERATOR EIGENVALUES (INCLUDING EXTRANEOUS ZERO VALUES): 
( 0.759912)+J( 0.000000) 
( 0.000000)+3( 0.000000) 


RESIDUES AT THE POLES: 


POLES EMS I DUES 
REAL (À) IMAG(B) 


( 0.711854)*J( 0.564795)  ( 9.841431) EXP(A*T)*COS(B*T) 
( 0.711854)«J(-0.564795) | (-0.837409) EXP(A*T)*SIN(B*T) 


..TRANFUNC COMPLETED........ FILES FOR OPTGRAPH GENERATED.. 
DO YOU WISH TO GENERATE AN "OPTMATD" DATA FILE TO SAVE THE: 
{A} (SYSTEM), (B) (CONTROL), 

{Н} (OBSERVABLES), (G) KMEASUREMENT- FEED FWD. DIST.), 
(GAMD) (NOISE), (F) (CONTROL GAIN), 

(K) (FILTER GAIN), {Q} (OUTPUT COST), 

{К} (CONTROL COST), (V1) (PROCESS NOISE INTENSITY (PSD)) 


(V2) (MEASUREMENT NOISE INTENSITY (PSD)), 
(S) (DISCRETE CROSS PRODUCT WEIGHTING), 
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MATRICES FOR REENTRY TO THE OPTSYSX OR ORACLS PROGRAMS? 
(X OREN) 


=» = че еле мо е  . mí тыз сы. мо е е —= е се е е се е е се ем ме се —_ се со е се е се - се е и е е е е се се е се е е оо е се см е се ‘U s. 


ч еә =з == =з =з» мо мес ее ее се е сем --. ш= е е е см е е е см оо е е ое е е е е е е =з =з =з т=ф мо е =з =з -. =» ч» оо е “ә мо мош тыз -- < 
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